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OVERVIEW  AND  SUMMARY 


This  is  a  system  concept  study  of  how  the  VSC  Seismic 
system  could  operate  as  a  test  bed  for  developing  an  optimum 
automated  world-wide  seismic  network.  Many  problems  are  an¬ 
ticipated  in  developing  such  a  network.  These  are  discussed  in 
a  theoretical  way  and  algorithms  are  suggested  which  could  be 
tested  by  the  VSC  system.  These  would  lead  to  more  optimum 
performance  of  an  automated  seismic  network. 

The  VSC  Seismic  system  is  expected  to  realistically  emu¬ 
late  the  operation  of  a  seismic  surveillance  network.  This 
can  be  done  by  executing  a  set  of  candidate  algorithms  to  per¬ 
form  all  of  the  automatic  functions  which  transform  continuous 
seismic  sensor  data  into  desired  seismic  event  information. 
Operation  of  the  VSC  Seismic  System  is  expected  to  provide  a 
test  bed  for  objectively  evaluating  various  competing  automatic 
algorithms  which  detect,  locate,  and  identify  explosions.  It 
is  desirable  to  impelment  a  baselined  VSC  Seismic  System  with 
sufficient  flexibility  to  test  algorithms  needed  for  optimum 
operation  of  a  seismic  surveillance  network.  For  that  pur¬ 
pose,  users  can  apply  automatic  algorithms  to  process  standard 
seismic  data  files  and  gauge  objectively  the  impact  of  a  pro¬ 
posed  algorithm  on  the  operation  of  a  seism;. c  network.  By- 
testing  these  algorithms  under  such  realistic  operation  con¬ 
ditions,  these  algorithms  under  such  realistic  operating  con¬ 
ditions,  the  best  performing  algorithms  can  then  be  selected 
for  implementation  in  data  centers  which  monitor  test  ban 
treaties  . 
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lts  soa^»  an->'  seismic  surveillance  system  transforms 
continuous  seismic  sensor  data  to  a  set  of  located  and  timed 
events  classified  as  earthquakes  or  explosions,  with  sufficient 
information  to  support  such  a  conclusion.  An  automated  system 
ror  achieving  this  goal  can  be  represented  topologically  as  a 
=eries  of  function  processes.  This  is  a  general  representation 
of  che  s>  stem,  and  it  is  independent  of  the  geographical  dis¬ 
tribution  of  the  sensors  and  location  of  the  data  processing 
elements  of  the  system.  Whether  or  not  the  data  processing  is 
centralized  or  distributed  throughout  the  network,  a  sequence 
of  functions  must  be  carried  out  to  reduce  raw  sensor  data  to 
the  desired  event  information. 


At  the  front-end  of  the  system,  where  continuous  sensor 
data  are  accessed,  an  automatic  detector  generates  a  set  of 
apparent  arrival  times  of  signals.  A  location/association 
process  transforms  the  asynchronous  stream  of  signal  detections 
into  an  event  stream.  This  provides  a  preliminary  description 
of  the  location,  origin  time,  and  characteristic  of  each  seismic 
event.  These  are  the  presumed  events  which  are  generated  by 
operating  the  system.  Given  the  event  locations,  an  automatic 
signal  editing  process  transforms  stored  continuous  sensor  data 
into  tentative  seismic  event  records.  These  records  contain 
signals  associated  with  an  event  or  background  noise  at  the  ex¬ 
pected  arrival  of  the  phase  of  interest.  At  this  stage,  records 
rrom  arrays,  three- -component  sensors,  or  s ingle  -  sensors  are 
reduced  to  single- - channel  records  of  seismic  phases.  These 
are  composed  or  compressional ,  shear,  or  surface  wave  phases. 
Next,  a  detector,  optimized  to  accurately  time  each  seismic 
phase,  activates  processes  to  measure  signals  or  noise,  and  ex¬ 
tracts  a  short  compressed  edit  of  the  detected  seismic  phases. 
‘aese  signal  measurement  data  and  compressed  seismic  phase 
edits  are  subsequently  input  to  an  interactive  partition  of 
tr.e  V  SC  seismic  System.  This  performs  functions  to  generate 
event  discriminants,  refine  the  location,  and  classify  the 

event  e.g.,  as  an  earthquake,  an  explosition,  or  an  unknown 
event  zvne)  . 
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ihe  concept  of  a  sequence  of  linear  programmed  function 
processes,  each  of  which  independently . access  and  updates  data 
rtxes,  is  illustrated  by  Figure  5-1.  The  decision  functions 
which  control  the  flow  of  information  are  shown  in  diamond  shape 
coxes.  The  functions  to  be  performed  are  shown  in  rectangular 
coxes.  Major  data  files  stored  by  the  system  are  also  shown. 
Major  data  files  stored  by  the  system  are  also  shown.  An  as¬ 
terisk  is  placed  alongside  of  those  functional  processes  covered 
oy  our  Automatic  Seismic  Signal  Processing  Research. 

Inspection  of  Figure  S-l  makes  clear  that  certain  tradeoffs 
are  involved  in  designing  the  VSC  Seismic  System.  Clearly 
Figure  S-l  does  not  constitute  a  design  of  the  system,  but 
merelv  shows  the  basic  functions  to  be  performed  by  any  seismic 
surveillance  system.  It  also  shows  the  sequential  nature  of 
these  function  processes  and  their  interaction  with  data  file 
structures  generated  by  the  system.  For  example,  the  continu¬ 
ous  waveform  file  might  be  either  centraliced  or  distributed 
in  some  manner  throughout  the  network.  In  any  case,  the  signal 
editor  will  have  a  need  to  access  whatever  data  are  available 
in  order  to  make  signal  measurements.  The  design  of  the  VSC 
beismic  System  obviously  will  place  limits  on  performance: 

•  Communications  affect  the  amount  of  data  which  can 
be  assessed,  the  time  delay,  and  the  amount  of 
storage  required  for  data  accessed. 

•  The  complexity  and  reliability  of  the  automatic 
algorithms  also  affect  the  amount  of  raw  data 
storage  needed  and  time  delays  inherent  in  pro¬ 
cessing  the  data  . 

•  The  effectiveness  of  front-end  detection  and  loca¬ 
tion  afreets  the  amount  of  interactive  processing 
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required  for  quality  control.  This  indirectly  af¬ 
reets  the  amount  of  raw  data  which  needs  to  be 
stored  and  time  delayed. 

The  setting  of  standards  on  record  sice,  number  of 
phases  to  be  retrieved,  and  digital  sample  rates 
afreet  the  reliability  of  the  signal  editing  pro¬ 
cess.  It  also  affects  the  amount  of  storage  and 
computer  power  needed  for  automatic  processing. 

i  ..e  maximum  reduction  of  data  by  automatic  process¬ 
ing  reduces  the  modest  proportions  the  computer 
power  needed  for  interactive  processing. 


An  effective  phase  detector/ timer  and  compressed  edit 
process  greatly  reduces  the  amount  of  data  accessed 
bv  interactive  processes.  This  greatly  reduces  the 
amount  of  raw  data  storage  needed  and  the  time  delay 
of  the  system  in  reporting  results. 


Given  this  system  context  of  functional  requirements,  we 
will  design  an  automatic  seismic  signal  processing  package  which 
"ill  incorporate  some  existing,  proven  algorithms  for  signal 
editing,  extraction,  and  measurement  and  suggest  alternative 
algorithms  for  future  testing.  Moreover,  we  have  also  redesigned 
'~,'e  ^nSer  detec  .or  to  correct  problems  which  have  caused  it  to 
nuss  a  large  percentage  of  signals  in  the  past.  These  algorithms 
"iH  De  designed  tor  modular,  independent  operation  within  the 
VSC  Seismic  System  context.  As  such,  they  will  be  designed  for 
optional  execution  by  users  of  the  VSC  Seismic  System,  or  thev 
can  3e  replaced  by  alternative  algorithms  for  comparative  test¬ 
ing.  We  anticipate  that  this  goal  is  compatible  with  the  con¬ 
cept  of  a  VSC  Seismic  System  as  being  a  test  bed  for  examining 
networx  operation.  It  is  also  anticipated  that  other  algorithms 

put  into  the  baseline  of  the  VSC  Seismic  System  will  be  similarly 
structured . 
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DESCRIPTION’  OF  THE  BASELINE  AUTOMATIC  SIGNAL  EDITOR 
CASE) 


xhe  oaseline  from  which  the  Automatic  Signal  Editor 
LASE;  was  derived  was  the  system  used  by  ENSCO,  Inc.,  in  the 
Event  Identii ication  Experiment.  Its  purpose  and  the  proce¬ 
dures  it  utilized  to  edit  seismic  signals  will  be  briefly  re¬ 
viewed  below.  The  Automatic  Signal  Editor,  as  we  define  it,  is 
cr.aL  part  of  a  seismic  surveillance  system  which  independently 
processes  and  derives  measurements  from  seismic  events  and  up¬ 
dates  files  of  such  measurements  for  use  in  the  interactive 
(analyst  dominated)  portion  of  the  seismic  processing  system. 


1.  Purpose  and  Generalized  Procedures  Associated  with 
the  ASE 


The  purpose  of  an  advanced  Automatic  Signal  Editor  is 
to  selectively  improve  preliminary  seismic  source  estimates, 
including  source  location,  origin  time,  depth,  magnitude, 
and  discrimination  parameters,  as  well  as  to  reduce  the  number 
of  unassociated  or  misassoc iat ed  seismic  phases  which  are  de¬ 
tected.  These  efforts  must  be  accomplished  prior  to  initiation 
ot  the  analyst  intensive  interactive  processing  phase  of  analysis 
to  insure  a  time  efficient  interactive  session.  In  concept 
The  ASE  could  be  employed  as  follows. 

After  an  event  has  been  tentatively  located  and  its 
origin  time  estimated  by  an  automatic  association  program,  the 
arrival  time  of  selected  seismic  phases  can  be  predicted  at 
each  available  station.  A  timing  tolerance  is  established 
tor  worst  case  location  errors,  allowing  enough  time  to  ob¬ 
serve  signals  trom  complex  source  regions  with  multiple  propa¬ 
gation  paths.  Sufficient  time  is  also  allowed  to  accurately 
determine  the  noise  state  prior  to  the  arrival  of  the  signal. 
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t  Given  a  set  o£  records  containing  signals  of  unknown 

precise  start  tine  and  dura t  ion ,  automatic  processors  are 
,  used  t0  U)  select  records  containing  useable  seismic  data, 

L-  condition  data  to  minimize  the  effect  of  malfunctions, 
ana  (5)  to  utilize  an  automatic  detector  algorithm  to  more 
precisel)  time  the  onset  of  the  seismic  signal  and  determine 
its  duration. 

Those  portions  of  the  seismic  record  which  contain  one 
or  more  possible  seismic  signals  from  the  event  are  extracted 
°m  record  and  filed  as  possible  signal  waveforms  of  the 
known  seismic  event.  These  siesmic  waveforms  can  be  accessed 
later  by  seismic  analysts  for  quality  assurance,  reprocessing, 
or  po s t - do t ec t io n  validation  of  their  presumed  signal  status. 

The  bottom  line  of  the  signal  editing  process  is  to  mea¬ 
sure  the  character  of  the  seismic  signal.  The  signal  measure¬ 
ments  are  filed  for  later  access  by  the  analyst.  They  can  be 
|  used  for  research  and  to  perform  the  functions  of  event  dis¬ 

crimination  and  characterization, 
i 

2.  Descripton  of  Seismic  Event  Editing  Procedures 

Continuous  seismic  waveform  data  are  accessed  to  extract 
segmented  records  containing  desired  signals.  These  are  short- 
period  regional  P  waves,  S  waves,  short-period  surface  waves 
and  long-period  surface  waves. 

An  automatic  detector  algorithm  is  applied  to  extract 
more  precise  time  and  measure  the  magnitude  of  signals  from 
known  events.  Cur  experience  with  Unger's  algorithm  in  the 
Event  Identification  Experiment  has  shown  that  this  procedure 
currently  severely  limits  the  effectiveness  of  the  signal  edit¬ 
ing  process.  As  a  result  there  is  a  serious  performance  gap 
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retween  the  manual  retrieval  of  seismic  signals  by  seismic 
analysts  and  that  of  automatic  detecors.  Our  goal  is  to  over¬ 
come  this  performance  gap  by  redesigning  the  automatic  de¬ 
tector  which  times  and  extracts  seismic  signals  from  long  seis¬ 
mic  records  . 


B.  DESCRIPTION  OF  UNGER'S  DETECTOR 
1 .  Theoretical  Background 

Tne  general  model  for  seismic  noise  consists  of  random- 
1>  modulated  envelope  and  phase  angle  functions.  The  noise 
statistics  underlying  the  modulation  are  presumed  to  be  sta¬ 
tionary  over  the  duration  of  the  seismic  record.  Therefore, 
noise  parameters  which  are  derived  from  noise  preceding  the 
signal  can  be  applied  over  the  entire  duration  of  the  record. 

The  general  model  for  seismic  signals  is  a  randomly 
modulated  envelope  and  phase  which  starts  at  the  Kth  point 
and  ends  at  the  CK*D)th  point.  With  respect  to  retrieving 
seismic  signals,  the  detection  problem  is  to  correctly  decide 
wnere  the  signal  starts,  i.e.,  on  the  Kth  point,  and  how  long 
it  lasts,  i.e.,  for  a  duration  of  D  seconds. 

Analytical  detectors  are  based  on  amplitude  and  phase 
angle  modulation  measurements  which  are  subject  to  error.  The 
rrequency  of  a  signal,  measured  as  the  time  derivative  of  the 
phase  angle,  is  particularly  subject  to  large  errors  due  to 
tne  large  envelope  fluctuations.  Considerable  effort  is  needed 
t0  obtain  more  accurate  and  precise  measurements  of  a  signal's 
frequency  so  that  the  frequency,  as  well  as  magnituae, can  be 
used  for  character  recognition  of  weak  signals. 

1.  Unger's  Criteria  for  Timing  Signals 
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Unger  greatly  simplified  the  theory  of  noise  by  consider¬ 
ing  noise  interterence  with  a  coincident  signal  as  a  constant 
level  of  power  subject  to  random  phase  modulation.  He  con¬ 
sidered  seismic  signals  as  a  fixed  amplitude  level  combined 
with  a  deterministic  phase  angle  versus  time  relationship.  For 
example,  a  turned-on  cosine  function  can  be  described  as  a  fixed 
amplitude  of  one  combined  with  a  linear  phase  versus  time  re¬ 
lationship  over  the  duration  of  the  signal. 

Unger  derived  a  relationship  for  determining  the  prob¬ 
ability  that  the  envelope  of  a  signal  combined  with  noise  in¬ 
terference  exceeds  the  fixed  noise  envelope  level;  and  that 
phase  angle  prediction  errors  are  less  than  tt/2.  He  observed 
that  detectors  based  on  this  phase  angle  criteria  can  certainly 
detect  signals  at  a  level  6  dB  lower  than  by  the  envelope  cri¬ 
teria.  However,  since  he  could  not  derive  a  deterministic  model 
for  the  time  variation  of  the  phase  angle,  he  could  not  success¬ 
fully  utilize  the  phase  angle  relationship. 

Unger  utilized  the  envelope  probability  relationship  to 
time  the  onset  of  seismic  signals.  He  computed  the  fraction 
of  times  a  possible  signal  exceeds  the  maximum  observed  noise 
and  applied  a  threshold  to  the  detection  statistic. 

3.  Design  of  Unger's  Analytic  Detector 

Unger  determined  the  peak  level  of  noise  preceding  the 
portion  of  the  record  containing  possible  seismic  signals. 

He  applied  an  envelope  probability  threshold  of  0.3  to  a  for¬ 
ward  looking  time  gate  of  duration  4  seconds.  That  is,  if 
he  found  that  a  maximum  fraction  of  envelope  values  exceeded 
tne  maximum  values  observed  of  noise,  he  would  declare  a  pos¬ 
sible  signal.  Then  he  searched  for  the  first  signal  peak  at 
least  2  to  3  dB  above  the  maximum  noise  to  confirm  the  signal 
detection  and  time  its  onset.  Finally  he  would  back  up  3/4 
of  a  cycle  to  precisely  time  the  signal  onset. 
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Poit-Mortem  Evaluation  of  Unger's  Algorithm 
Applied  as  an  Automatic  Signal  Editor 


In  :he  Event  Identification  Experiment,  Unger's  algorithm 
•vas  used  to  time  the  occurrence  of  signals  from  known  events 
on  long  time  records.  It  controlled  the  decision  as  to  what 
portion  of  a  short-period  record  should  be  extracted,  measured, 
and  filed.  The  detector  was  part  of  a  larger  automated  Short- 
Period  Earthquake  Editor  (SPEED)  package  which  accessed  records, 
corrected  for  system  response  and  preconditioned  the  long  records 
01  data  prior  to  extracting  and  measuring  signal  waveforms, 
correspondingly  in  the  Event  Identification  Experiment,  a  Long- 
Period  Earthquake  Editor  (LPEED)  performed  the  same  function 
for  long-period  surface  wave  data. 


Two  aspects  of  the  performance  of  Unger's  detector  were 
excellent.  The  timing  precision  of  the  detected  signals  was 
estimated  to  be  approximately  +  0.1  seconds  (about  the  level 
which  would  be  expected  from  seismic  analysts).  Also,  the 
false  alarm  fraction  of  extracted  signals  was  low  (about  5%)  . 
In  this  respect  it  exceeded  acceptable  performance  (a  false 
alarm  rate  of  10%)  .  Lowering  the  threshold  to  detect  more 
signals  was  attempted  but  did  not  significantly  change  the 
detection  capability  of  the  algorithm. 

For  acceptable  performance,  an  automatic  signal  editor 
should  be  able  to  detect  most  of  the  signals  visible  to  an 
analyst.  An  acceptable  level  would  be  a  detection  probability 
Pd  1  Unger's  automatic  detector  applied  to  signals  edited 

by  seismic  analysts  was  not  capable  of  achieving  this  level  of 
performance  at  any  false  alarm  rate.  Unger's  algorithm  was 
unable  to  detect  and  time  any  signals  for  approximately  40% 
of  the  events  edited.  This  indicates  a  detection  probability, 
Pd  <  °-°-  rhe  detector,  in  its  present  form,  has  a  missed 
signal  problem  and  needs  to  be  revised  to  attain  an  acceptable 
j.evel  of  signal  extraction  performance. 
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another  requirement  tor  acceptable  performance  is  to 
'.--icate  signals  extracted  by  the  automatic  detector.  Net- 
•vcr.\  \  a. i cation  couic  be  achieved  by  extracting  signals  with 
extreme  magnitude  and  arrival  time  errors  .  In  its  present 
-o.m,  Spoil)  does  not  carry  out  this  network  validation  f unc - 
Consequently,  it  occass  ionally  detected,  timed,  and. 
measured  mixed  event  signals,  which  although  obviously 
MsiD.e  to  an  analyst  were  not  valid  because  of  apparently 
extreme  magnitude  and  travel  time  anomalies.  We  estimated 

t. cat  o-j  to  12%  or  the  long  records  contain  mixed  event  sig¬ 
nals.  Clearly  this  type  of  error  must  be  reduced  to  negli- 
3-b_e  proportion  since  it  can  result  in  erroneous  event 

u. iaracterizat ion ,  especially  for  smaller  events. 


Other  major  raults  were  observed  in  our  use  of  Unger's 
detection  algorithm  in  SPEED  for  the  Event  Identification 
Experiment.  The  algorithm  missed  about  50%  of  the  signals 
detected  by  analysts,  which  were  visible  to  an  analyst  but 
cic  not  exceed  the  maximum  observed  noise  level.  The  algo- 
rit.tm  missed  another  10  to  151  of  signals  which  were  obvious 
-z  an  analyst  and  clearly  well  above  the  maximum  observable 
s -  I'-vel.  Our  post-mortem  assessment  of  errors  indicated 
22\era_  causes.  Unger's  method  of  determining  the  maximum 
noise  level  was  not  robust.  Occasional  spikes  or  glitches 
ir.  one  noise  caused  his  algorithm  to  set  an  unreasonably 
..igh  maximum  noise  level  shich  could  not  be  exceeded  even  by 
'  a*“-’  visible  signals.  Unger's  model  showed  a  tendencv  to 

miss  some  clearly  visible  impulsive  signals  of  short  dura¬ 
tion.  The  cause  of  this  effect  was  the  use  of  a  4  second 
signal  duration  for  estimating  the  probability  that  maximum 
wws  e.  .  ec  ac.se  levei  is  exceedec.  Impulsive  signals  of  less 
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1.3  seconds  duration  would  fail  to  be  sensed  by 
the  threshold  criteria  applied  to  the  estimated  probability 
or  exceeding  maximum  noise.  Impulsive  signals  are  not  un¬ 
common  and  it  is  thus  essential  that  the  presumed  duration 
or  signals  cover  a  much  wider  range,  e.g.,  from  0.3  seconds 
to  10  seconds.  Unger's  detector  missed  some  complex  signals 
tor  mucn  the  same  reason.  These  extend  over  a  time  span 
greater  tnan  4  seconds  and  the  envelope  of  these  signals 
are  nignly  variable  and  skewed  toward  low  level  fluctuations 
.n  tnose  cases  pulses  above  the  maximum  noise  level  are  fre¬ 
quently  interspersed  with  dead  spots  below  that  level.  This 
results  m  a  failure  of  the  probability  test  to  stop  the 
cetector  at  the  beginning  of  such  complex  signals.  These 
are  commonly  missed  or  detected  with  very  late  start  times. 
Unger's  algorithm,  as  presently  implemented  in  SPEED,  warms 
up  on  one-minute  of  noise  preceding  the  window  containing 
signals  and  continuously  updates  noise  statistics  until  a 
signal  is  detected.  This  updating  procedure  reduces  the 
robustness  of  the  algorithm,  especially  in  the  case  of  de¬ 
tecting  and  timing  emergent  events.  In  that  case,  a  signal 
wnicn  is  initially  less  than  the  maximum  noise  level  rises 
gracualiy  above  that  level  but  not  fast  enough  to  pass  the 
probability  threshold  test.  Instead,  the  maximum  observed 
r.cise  level  is  gradually  raised  and  effectively  shuts  off 
tr.e  detector  so  that  it  misses  emergent  signals.  Another 
rault  which  caused  major  errors  in  subsequent  event  identi- 
-ication  was  the  occasional  erroneous  extraction  of  the 
wrong  event.  Since  event-signals  are  accessed  with  four 
minute  time  windows  and  events  can  be  expected  at  the  rate 
or  about  one  per  hour,  there  is  a  small  but  significant 
cnance  ^arcut  o  to  10%)  that  more  than  one  event's  signals 
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sre  accessed  by  such  a  long  record.  Since  SPEED  has  no  post- 
dececticn  validation  tests  of  extracted  signals,  this  effect 
results  in  an  occasional  serious  editing  problem,  especially 
for  small  events  which  are  interfered  with  by  larger  events . 
The  final  fault  of  SPEED  is  that  it  is  not  fully  automatic. 
Records  were  visually  scanned  by  us  to  remove  obvious  mal-  • 
functioning  stations.  In  a  fully  automated  signal  editor, 
this  quality  control  should  be  built  into  the  signal  editor. 

As  a  result  of  our  post-mortem  evaluation  of  SPEED,  we 
are  proposing  major  revisions  to  obtain  a  fully  automated 
signal  editor  which  has  the  potential  of  performing  at  a 
level  comparable  with  that  of  a  seismic  analyst.  In  order 
to  time  and  measure  weak  signals  at  levels  below  the  most 
probable  occurring  noise  level,  we  propose  applying  a  multi¬ 
variate  analytic  detector  capable  of  sensing  the  character 
of  such  signals,  i.e.,  frequency,  bandwidth,  duration,  etc. 
as  well  as  sensing  positive  power  fluctuations.  The  multi¬ 
variate  analytic  detector  should  be  designed  to  detect,  time, 
and  determine  the  duration  of  multiple  arrivals  at  each  sta¬ 
tion  to  be  sorted  out  and  interpreted  by  a  pos t- detection 
network  validation  strategy.  Signal  and  noise  parameters 
sr.culd  be  measured  by  robust  ordered  statistical  analysis 
wnich  is  insensitive  to  occasional  spikes  and  other  inter¬ 
mittent  malfunctions .  No  prior  assumption  should  be  made 
of  the  signal  duration.  This  parameter  is  highly  signal 
dependent  and  should  be  determined  effectively  and  efficient¬ 
ly  from  measurements  of  data.  There  is  no  need  to  continu¬ 
ously  update  noise  statistics  as  is  dene  for  front-end  de¬ 
tectors;  noise  parameters  should  be  assessed  from  noise 
preceding  the  signal  and  assumed  to  apply  over  the  subsequent 
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curat  ion  or  :r.e  record.  Network  validation  procedures  need 
to  be  applied  to  assure  that  any  signal  extracted  is  rea- 
scr.aoly  consistent  with  what  is  expected  from  the  known 
event.  Automated  quality  control  will  be  performed  to  re¬ 
move  obvious  malrunctiomng  stations  and  to  condition  mar- 
gma.  data  subject  to  intermittent  malfunctions.  Finally,  . 
a  reasonably  acceptable  automatic  signal  editor  should  ef- 
rectively  trade-orr  false  alarms  for  improved  detection  of 
weak  signals  from  small  events. 


C.  ANALYSIS  OF  AN  ADVANCED  ANALYTIC  DETECTOR 

1.  Generalized  Analytic  Signal  and  Noise  Model 

ihe  ambient  noise  model  will  be  generalized  to  include 
the  effects  of  system  malfunctions.  Based  on  Unger's  [1978) 
study  of  noise  envelope  fluctuations,  normal  noise  envelope 
r luctuations  are  presumed  to  be  consistent  with  or  close  to 
a  Gaussian  noise  model.  System  malfunctions  will  appear  as 
severe  departures  of  extreme  values  from  such  a  model.  This 
generalized  noise  model  provides  a  basis  for  building  auto¬ 
matic  quality  control  into  the  automatic  signal  editor. 


The  generalized  signal  model  also  will  have  a  system 
malrunction  component  to  avoid  contamination  of  automatic 
signal  measurements  by  occasional  intermittent  system  mal- 
r unctions .  The  signal  model  is  generalized  beyond  that  of 
Unger's  (which  presumes  a  single  signal  arrival)  to  include 
multiple  signal  arrivals.  These  will  include  multiple  sir- 
ra_s  rrom  complex  events  (i.e.,  multiple  explosions  or 
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earthquakes.',  secondary  propagation  phases,  regional  phases, 
v-..e  pnases,  nixet  event  signals  and  local  events,  and 
.  wC —  e.  sua^ered  phases  unique  to  a  particular  seismic  sta- 
.ne  system  noise  model  will  include  effects  caused  by 
sr_.<es  =.nd  gli-cnes,  dead  spots  and  clipDing. 


Multivariate  Amplitude  and  Phase  Statist! 


cs 


We  iefine»  on  exacting  physical  basis,  a  set  of 
2  v-a. w  mul -i v  ariate  signal  amplitude  and  frequency  measure¬ 
ments  .  These  are  designed  to  achieve  band-limited  whitening 
c~  ^e-smic  signals  or  any  type  which  can  be  reasonably  expect¬ 
ed.  Spectral  moment  operators  are  applied  for  this  purpose. 
Measurement  or  variates  related  to  the  bandwidth  of  observed 
signals  and  measurement  of  the  dominant  frequency  of  observed 
signals  are  used  for  feed-back  control  of  post-bandpass  fil- 

'■w*s  improve  extraction  of  band-limited  pre-whitened  seis¬ 
mic  signals. 

Basic  Ground  Motion  Measurements 
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.-.n  important  r actor  in  implementing  Unger's  theory  of 
5e  -electors  is  the  error  and  stability  problem  of  com- 
.n5  accurate  phase  and  rrequency  measurements.  The  prob¬ 
es  caused  by  the  interdependence  between  random  frequency 
a-iu-li tude  modulations  or  seismic  noise  and  signals.  To 
■*e  tms  problem  we  will  transform  the  data  to  an  asynchron- 
time-ordered  set  of  envelope  peaks.  On  that  basis,  our 
.ysis  indicates  that  we  can  condition  the  data  to  obtain 
:  more  accurate  and  precise  measurements  of  the  magnitude 
; round  motion,  frequency,  arrival  time,  and  the  pulse 
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width  associated  with.  each,  envelope  peak.  In  order  to  achieve 
o and- United  pre-whitening  of  any  conceivable  type  of  signal, 
w e  will  apply  five  spectral  moment  (^difference  operators)  for 
that  purpose.  In  this  way  any  plausible  seismic  source  is 
automatically  searched  for  whitened  bands  where  the  source 
can  be  accurately  measured.  For  example,  if  the  frequency, 
falls  off  as  f--  above  the  corner  freauencv  then  a  second 


derivative  operator  with  spectrum  f 
spectrum  above  the  corner  frequency 


.  ? 


whitens  the  source 


Interpretation  of  3asic  Ground  Motion  Measurements 
of  Noise;  Data  Conditioning 


From  the  analysis  of  Rice  [1954),  we  can  detect  the 
presence  of  strong  almost  periodic  noise  interference  as 
significant  departures  from  expected  Gaussian  noise  statis¬ 
tics.  A  technique  was  derived  for  demodulating  such  noise 
to  separate  the  desired  broadband  data  from  the  interfering 
nearly  periodic  noise.  Where  such  problems  are  encountered 
we  will  demodulate  the  desired  data  from  such  interfering 
components . 


Based  on  a  relationship  observed  by  Unger  relating  the 
frequency  of  seismic  envelope  modulations  to  the  standard 
deviation  of  phase  measurements,  we  expect  to  see  a  large 
dependence  between  the  magnitude  fluctuations  of  seismic 
"rise  or  signals  and  the  frequency  or  spectral  bandwidth  of 
such  envelope  peaks.  A  theory  for  correcting  noise  for  ob¬ 
served  magnitude  and  frequency  fluctuations  w as  developed. 
The  result  of  applying  such  corrections  to  noise  data  is 
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ignif i cantly  reduce  the  variance  of  observed 
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requency  fluctuations  of  noise.  This,  combined 
lemoval  of  periodic  noise  components  is  expected  to 
s  iqmr  icantly  enhance  the  detection  of  weak  signals  by  inde- 
.e.iw.nt.}  a.-pl>  ing  magnitude  and  frequency  measurements. 

ihe  application  of  such  corrected  independent  detection 
statistics  will  be  applied  to  the  five  moment  operators  de- 
3-t.ned  ..or  band-limited  whitening  of  seismic  signals.  This 

’"in  result  in  the  definition  of  a  ten  component  multivariate 
measurement  of  seismic  envelope  peaks.  In  general,  and  es¬ 
pecially  for  weak  signals  it  is  not  expected  that  all  of  the 
components  will  be  detected.  Those  detected,  however,  will 
oe  measured.  Those  not  detected  will  be  classified  as  noise 
which  will  be  measured  as  control  for  subsequent  network 

analysis  of  signal  measurements  (i.e.,  to  obtain  unbiased  es¬ 
timation  of  magnitude) . 


3.  Detection  of  Apparent  Signals 

A  theoretical  treatment  is  given  of  the  problem  of  re¬ 
ducing  basic  ground  motion  measurements  to  statistically  in¬ 
dependent  magnitude  and  frequency  fluctuations.  The  input 
aata  to  such  a  process  are  observed  magnitude  fluctuations 
corrected  to  remove  any  frequency  dependence  and  observed 
frequency  fluctuations  corrected  to  remove  any  magnitude  de¬ 
pendence.  A  linear  rotational  transformation  operator  Q  is 
applied  to  transform  1  statistics  of  corrected  magnitude 
fluctuations  and  :  statistics  of  corrected  frequency  fluc¬ 
tuations  to  independent  unit  variance  measures  of  magnitude 
fluctuations  of  noise  and  frequency  fluctuation  of  noise. 
Significant  deviations  of  signals  from  these  di stributions 
provide  an  independent  basis  for  detecting  a  signal  based  on 

i " “  ClianSe  °f  frequency  or  its  change  of  magnitude  from 
values  expected  from  noise. 
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?uiie  width  measurements  of  noise  are  related  to  the  measured 

<-i.te  interval  between  noise  envelope  peaks.  Rice  (1954) 
s..ouei*  tiiai  this  statistic  can  be  used  to  estimate  the  bandwidth 
of  an  ideally  filtered  Gaussian  process.  By  performing  this 
calibration  we  can  use  our  basic  pulse  width  measurement  of 
envelope  peaks  to  estimate  the  bandwidth  of  a  Gaussian  process 
whitened  by  the  spectral  moment  operators.  This  is  extremely 
vaiuabie  information  because  the  observed  dominant  frequency  ' 
ana  tne  associated  bandwidth  provide  a  basis  for  feedback  con¬ 
trol  or  a  bandpass  filter  optimally  designed  to  extract  the 
band-nmited  whitened  Gaussian  process.  This  will  be  applied 
to  marginal  detections  of  weak  signals  to  confirm  their  de¬ 
tection  status  and  to  extract  those  signals  with  minimal  seis¬ 
mic  noise  interference. 

b.  Detecting  Complex  Signals 

The  problem  of  detecting  signals  as  a  single  significant¬ 
ly  unusual  envelope  peak  is  generalized  to  that  of  detecting 
strings  of  such  outlier  envelope  peaks.  The  latter  envelope 
peaks  would  be  expected  to  be  replicas  of  the  same  signal  pro¬ 
cess  which  would  be  expected  to  yield  the  same  magnitude  and 
r r equency  fluctuations.  Thus  by  ooserving  repeated  anomalous 
peaks  which  vary  little  from  each  other,  we  can  apply  the  like¬ 
lihood  aetection  criteria  proposed  by  Unger  (1978)  which  should 

ennance  the  detection  of  very  weak  but  persistent  complex  sig¬ 
nals  . 


ihe  cetection  of  signals  as  strings  of  envelope  peaks 
provides  additional  event  information.  The  start  time  of 
tne  signal  is  taken  as  the  measured  start  time  of  the  first 
peak  in  the  cluster  of  peaks.  The  end  time  is  taken  as  the 
ena  time  of  the  last  envelope  peak  of  the  signal  process. 
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Measurements  of  complex  signals  will  be 
orcerea  statistical  analysis  of  envelope  peak 
tuations  and  frequency  fluctuations  within  the 
containing  the  clustered  envelope  peaks. 


accomplished  by 
magnitude  fluc- 
time  window 


Network  Strategy  for  Retriving  Signals 
small  events 


from 


■\e  expect  to  be  faced  with  the  problem  of  identifying 
small  events  which  are  masked  by  seismic  noise.  By  applying 
taxed  constant  false  alarm  rate  thresholds  to  constrain  false 
alarms,  we  may  keep  such  false  alarms  to  an  acceptably  low 
level,  but  at  the  same  time  make  it  nearly  impossible  to  ex- 
u*act  *eak  signals.  In  the  case  of  low  magnitude  signals, 
we  can  employ  a  station  variable  false  alarm  rate  strategy 
which  raises  the  expected  ratio  of  retrieved  signals  to  re¬ 
trieved  false  alarms  to  an  acceptable  level.  On  the  presump¬ 
tion  that  an  identification  is  required  of  all  events  analyzed, 
it  is  better  to  have  at  least  one  or  a  few  signals  combined 
with  some  false  alarms  than  to  have  only  noise  measurements 
to  identify  the  signal.  This  strategy  is  made  somewhat  more 
sensible  by  the  fact  that  network  validation  procedures  will 
06  aDlS  t0  eliminate  at  least  some  of  the  false  alarms  and 
tr.at  network  magnitude  determinations  will  compensate  to  some 
extent  for  noisy  determinations  of  the  event  parameters. 

D-  SUMMARY  AND  RECOMMENDATIONS 


.  A  conci3e  description  will  be  given  for  the  design  of  a 
multivariate  analytic  detector  used  to  extract  short-period 
seismic  signals.  From  this  we  will  be  able  to  gauge  the  mini¬ 
mum  errort  required  to  implement  this  srategv  for  retrieving 
signal ^information.  From  this  we  expect  to  develop  the  con¬ 
cept  of  a  baseline  implementation  of  a  multivariate  analytic 
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-etector  which  is  sufficiently  advanced  to  achieve  our  goal  of 
editing  all  possible  weak  signals  from  small  events.  The  out¬ 
put  or  an  automatic  seismic  signal  editor  is  to  provide  the 
maximum  amount  of  additional  data  from  known  seismic  events 

to  improve  location  of  the  event  and  to  improve  identification 
or  the  source  . 


One  or  the  functions  performed  by  the  automatic  signal 
eaitor  is  to  extract  all  of  the  waveforms  which  can  be  possi- 
o  >  interpreted  as  signals  from  the  event.  These  extractions 
include  long-period  surface  waves,  short-period  teleseismic 
-  waves  and  secondary  phases,  and  regional  phases.  This  in- 
rormation  will  be  automatically  inserted  into  signal  measure¬ 
ment  riles  which  can  be  accessed  for  interactive  seismic  pro- 
cessing  by  seismic  analysts. 

Another  function  of  the  automatic  signal  editor  is  to 
reduce  basic  ground  motion  measurements  to  estimates  of  the 
event  magnitude,  the  dominant  frequency  and  frequency  band, 
an  the  complexity  of  an  event  and  its  coda.  Some  advanced 
applications  of  the  signal  measurements  will  be  to  derive  source 
aiscr lmmatnts ,  application  of  clustering  theory  to  identify 
anomalous  events,  and  source  region  calibration  to  determine 
precise  magnitudes  of  normal  earthquake  events. 
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SECTION  I 

DESCRIPTION  OF  THE  BASELINE  FOR  DESIGNING  AN 
AUTOMATIC  SIGNAL  EDITOR  CASE) 


Our  background  of  performing  the  Automatic  Signal  Editing 
' "  -action  was  the  Event  Identification  Experiment,  Sax, 
e-  al.  (19/9).  There  we  transformed  continuous  waveform  data 
c‘  .ong  records  into  compressed  extractions  of  seismic  signals 
and  into  a  set  of  seismic  measurements  to  be  used  for  event 
icentif ication .  The  functions  which  were  performed  completely 
under  automatic  control  were  to: 


Access  seismic  records  containing  the  desired  seis 

signals  from  events  of  known  location  and  mag¬ 
nitude  . 

Utilize  an  automatic  detector  to  detect  and 
accurately  time  the  signals. 

Generate  files  of  compressed  seismic  signal  wave¬ 
forms  . 


•  Generate  files  of  signal  measurements. 

The  files  generated  by  the  AGE  were  then  accessed  by  an  inter- 
ac-*1^  seismic  processor  which  transformed  the  signal  measure¬ 
ments  to  a  set  of  event  measurements  characterizing  the  source 
-he  seismic  event.  From  the  event  measurements,  we  gen¬ 
erate!  discriminants  and  performed  adaptive  statistical  anal¬ 
yses  to  identify  obvious  earthquakes  and  to  identify  anomalous 
seismic  events  such  as  explosions. 


A. 


PURPOSE  AND  GENERALIZED  PROCEDURES  ASSOCIATED  WITH  THE  ASE 


Some  obvious  benefits  of  speed,  efficiency,  and  objec¬ 
tivity  are  obtained  by  automatically  reducing  the  flow  of 
continuous  waverorm  data  to  the  essential  data  required  to 
pertorm  event  identification.  Why  this  rapid  reduction  of 
cata  is  a  userul  goal  can  best  be  illustrated  by  the  follow¬ 
ing  scenario  of  world-wide  seismic  data  collection. 

•  Events  of  m^>4  are  expected  at  the  average  rate  of 
one  per  hour. 

•  A  network  employed  to  measure  the  events  includes 

10  short-period  (SP)  vertical  component  single¬ 
sensors 

S  nineteen  element,  SP,  vertical  component  ar¬ 
rays 

10  SP,  three-component,  regional  sensors 
25  long-period  CLP),  three-component  sensors. 

•  The  input  data  flow  of  this  network  is 

10  million  words  per  event  of  SP  data 
90  thousand  words  per  event  of  LP  data. 


Since  most  of  this  data  collected  by  the  network  is  seis¬ 
mic  noise,  considerable  reduction  of  data  is  possible  by  sav¬ 
ing  only  data  containing  visible  signals.  To  estimate  the 
scale  of  this  reduction  consider  that  an  event  would  typically 
consist  ot  12  detected  short-period  signals  per  event  and  4 
cetected  long-period  signals  per  event.  Since  signal  durations 
are  expected  to  vary  between  one  and  ten  seconds,  the  com- 
pressec  edit  or  each  signal,  allowing  sufficient  space  for 
coca  and  noise  preceding  the  signal,  can  be  accomplished  with 


o.OO  were,  recores  or  aata  sampled  at  20  words  per  second. 

This  assumes  that  an  automatic  detector  can  reliably  detect 
and  time  the  seismic  signals. 

The  transformation  of  continuous  waveform  data  to  a  set 
or  compressed  signal  records  results  in  a  reduction  of  data 
from  1C  million  words  per  event  to  10  thousand  words  per 
event  or  short-period  signal  data.  A  further  reduction  of 
cata  or  10:1  is  possible  by  reducing  each  detected  signal 
wO  a  set  or  signal  measurements  such  as  the  arrival  time, 
measurements  of  magnitude,  complexity,  etc.  This  latter 
.eduction  of  data  to  one  thousand  words  per  event  is  a  dif- 
r icult  step  to  take,  exclusively,  because  of  the  problem  of 
agreeing  on  data  processing  standards  and  the  need  for  wave- 
rorm  data  to  maintain  historical  data  files.  At  some  future 
uime,  at  least  for  events  of  lower  priority,  one  could 
an cicipate  reducing  the  event  data  base  exclusively  to  sig¬ 
nal  measurements . 


By  reducing  continuous  waveform  data  to  compressed  sig¬ 
nal  records  and  signal  measurements,  it  will  then  become 
.easible  to  perform  analys t- interactive  processing  on  large 
dav.a  bases.  The'se  will  either  be  world-wide  or  from  targeted 
regions .  Our  experience  with  the  Event  Identification  Experi¬ 
ment  indicated  that  the  critical  element  controlling  the 
cuaiity  of  results  obtainable  by  this  automated  approach  is 
L.ne  pertormance  or  the  automatic  detector  used  to  detect  and 


time  seismic  signals.  Given 
bility  for  fast  and  reliable 


this  goal  of  developing  a  capa- 
event  identification,  we  will 


concentrate  on  developing  suitable  procedures 


and  algorithms 


-0.  extracting  and  measuring  signals  from  known  events. 
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*ns  v-^  ays  ten  should  provide  sufficient  noise  informa- 
-_w.i  -  -  om  non- detect  mg  stations  for  networking  algorithms 
sucn  as  those  for  estimating  unbiased  event  magnitudes. 
.-~.cur_>  RMS  noise  measurements  are  sufficient  to  provide  a 
.SiSonaole  basis  ror  estimating  the  relative  noise  at  dif- 
rerent  stations  of  the  network.  Such  hourly  RMS  noise 
figures  could  be  routinely  provided  as  event  header  informa¬ 
tion  with  only  50  words  per  event.  It  would  be  preferable 
also  provide  rrequency  dependent  noise  estimates .  This 
uwUld  be  done  on  a  daily  basis.  Such  network  noise  infor¬ 
mation  could  be  most  efficiently  provided  by  autocorrelation 
unction  estimates  of  seismic  noise;  by  64  lags  at  153  words 
per  event. 


Providing  this  important 
than  200  words  per  event.  It 
the  amount  of  storage  required 
seismic  signal  data. 


noise  information  would  be  less 
will  not  add  significantly  to 
in  files  containing  historical 


This  reduction  of  the  seismic  data  flow  to  sets  of  de¬ 
tected  signals,  signal  measurements  associated  with  the  event 
and  hourly  network  noise  update  information  provides  all  the 
information  needed  to  rapidly  identify  earthquakes  and  to 
main -ain  historical  seismic  event  data  files  required  for 
future  event  studies.  The  purpose  of  this  study  is  to  ex¬ 
plore  the  critical  procedures  and  algorithms  needed  for  this 
automated  approach  to  signal  editing  and  measurement  for 
ruture  application  to  the  VS C  system.  We  will  describe  some 
°f  the  e*xist--ng  algorithms  used  in  the  Event  Identification 
experiment  with  which  we  found  problems.  We  will  describe 
tnose  problems  along  with  the  modifications  needed  for 
satisfactory  operation  of  the  automated  VSC  system. 


1 

I 
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3.  DESCRIPTION  OF  SEISMIC  EVENT  EDITING  PROCEDURES 
'  USED  FOR  EVENT  IDENTIFICATION 

.he  automated  editor  carried  out  four  basic  analysis 
procedures.  These  are  to: 

•  Access  signal  time-windows 

•  extract  the  signal  waveform 

•  Measure  detected  signals 

•  Measure  noise  at  non-detecting  stations. 

i  1.  Accessing  Signal  Time  Windows 

a.  Short-period  signals 

I  Short-period  signals  are  accessed  from  the  computed  ar- 

i  rival  time  of  events  located  and  timed  by  the  VSC  system, 

j  The  initial  location,  depth,  and  timing  of  the  event  may  lead 

to  substantial  errors  in  timing  the  arrival  of  seismic  signals. 
The  rront-end  process  for  performing  this  function  is  shown 
schematically  on  Figure  I-l.  Less  than  the  desired  precision 
in  locating  the  rocus  of  events  stems  from  the  following 
factors. 

•  Multiple  transmission  of  signals  due  to  Comdex 
source,  path,  or  receiver  effects 

•  Large  uncertainty  of  correctly  associating  the  same 
phase  out  or  sets  or  multiple  arrivals  at  4  or  more 
stations 

•  Occasional  mixing  of  event-signals 


i 


•  Phase  association  dirticulties  due  to  magnitude 
variation  oetween  stations  and  multiple  transmiss ions 

•  Noise  masking  some  of  the  multiple  transmiss ions 

•  Dirr iculty  of  timing  emergent  and  complex  phases. 

As  a  result  of  these  uncertainties  and  difficulties  in  as¬ 
sociating  pnases  and  making  preliminary  focal  determinations; 
a. so,  as  a  result  of  the  present  state-of-the-art  of  automatic 
associate on/ location  process;  occasional  large  event  location 
and  timing  errors  are  to  be  expected.  Therefore,  in  order  to 
design  a  robust  ASE,  it  should  be  made  capable  of  correctly 
.etr-eving  event-signals,  even  under  such  adverse  circums tances . 
Prooably  at  least  several  percent  of  the  preliminary  locations 
will  be  in  error  by  as  much  as  10°  epicentral  distance.  Thus 
m  searching  for  signals  with  possibly  large  location  errors, 
tne  initially  selected  time-window  should  be  at  least  2  minutes 
m  duration.  Also,  at  least  one  minute  additional  is  needed 
to  assure  sufneient  noise  preceding  the  signal.  Another  one 
minute  interval  is  needed  for  coda  following  the  signal.  Thus, 
we  recommend  that  the  editor  search  a  4  minute  time  window  for 
each  desired  signal.  This  was  the  time  interval  used  by  us  in 
t.ie  Event  Identification  Experiment.  This  larger  window  pre¬ 
vented  most  of  the  serious  problems  our  automated  system  would 
nave  encountered  in  retrieving  signals  of  mis  located  events. 

On  tne  other  hand,  searching  this  larger  initial  time  window- 
increases  tne  chance  of  encountering  mixed  events.  These  odds 
arS  esti:nate-  to  oe  about  3$  to  6 $.  As  a  consequence  we  will 
require  more  sophisticated  application  of  post-detection  in¬ 
formation  to  correctly  select  out  the  valid  event- signals . 

0ur  ?er=°rmance  in  automatically  retrieving  signals  for  the 
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-vsnt  -centification  Experiment  suffered  from  the  lack  of 
;uch  post-detection  processing  to  validate  each  retrieved 
ignai.  This  situation  needs  to  be  improved  and  will  be 
Liscussed  later  under  recommended  modifications. 


b.  Long-period  signals 


Due  uO  die  long  duration  of  long-period  surface  wave 
signals,  we  saw  no  need  to  search  for  them  with  a  detector, 
tms  is  Decause  the  predicted  signal  time  interval  is  small 
compared  to  the  duration  of  such  signals.  Analytical  rela¬ 
tionships  were  derived  for  the  broad-region  dispersion  of 
surrace  waves  by  Unger  C1978)  .  These  are  shown  in  Figure  1-2 
for  the  Asian  continent  and  North  American  continent. 

Although  some  minor  difference  is  indicated,  a  slightly 
broader  envelope  can  be  used  to  retrieve  the  narrowband  passed 
surrace  wave  groups.  In  the  Event  Identification  Experiment, 
we  determined  a  positive  signal  detection  status  to  be  12  dB 
above  the  mean  noise  envelope  preceding  the  predicted  signal 
Scare  time.  We  round  this  procedure  to  be  satisfactory  in 
almo- 1  all  cases  . 


Regarding  the  retention  of  long-period  surface  wave  data 
m  permanent  read-only  files  for  advanced  waveform  analysis, 
we  round  that  the  short-period  limit  for  extracting  useful 
signal  information  was  12  seconds  period.  For  an  average 
spicentral  distance,  approximately  900  seconds  of  data 
are  required  to  sample  Love  waves  and  Rayleigh  waves.  A 


u-gicai  sampling  period  of  3  seconds  most 
v.ne  bandwidth  required  for  wave  groups  of 
or  -onger.  Thus,  two  records  of  300  words 


adequately 
12  seconds 
each  are  s 


covers 
period 
uff icient 


«n 


to  treats  remanent  read-only  files  for  retention  of  long- 
period  waveforms .  Broadband  and  narrowband  filtered  noise 
level  measurements  of  detected  signals  would  be  provided  as 
header  information  provided  with  signals  so  that  users  can 
determine  the  detection  status  of  signals.  Network  noise 
information  can  be  provided  by  the  procedure  described  for 
short-period  signals. 

Our  experience  with  the  Event  Identification  Experiment 
indicated  that  long-period  surface  wave  signals  could  be 
successfully  extracted  from  the  data  stream  by  deterministic 
timing  of  the  signals  using  Unger's  broad-region  dispersion 
relationship.  The  software,  LPEED,  used  for  that  purpose  can 
probably  be  transferred  with  only  minor  modifications  to  the 
VSC  system. 


By  contrast,  snort-period  signals  are  of  short  duration 
compared  to  the  much  wider  time  window  initially  accessed. 

It  is  therefore  necessary  to  apply  statistical  decision 
criteria  to  time  and  extract  these  short  waveforms  from  that 
much  longer  time  window  containing  noise.  This  problem  is 
caused  by  the  need  to  extract  weak  signals  of  short  but 
unknown  duration  from  ambient  noise  by  the  diverse  and  complex 
character  of  signals  to  be  extracted,  by  multi-phase  trans¬ 
mission  of  signals,  and  also  by  strong  interference  due  to 
receiver  scattering.  To  cope  with  these  problems,  we  applied 
an  automatic  detector  to  the  task  of  timing  and  extracting 
those  short-period  waveforms  possibly  associated  with  a  given 
event.  This  will  be  described  in  the  next  subsection. 


I.  Application  of  Automatic  Detectors  to  Extract  the 
Short -Period  Waveforms  of  Known  Seismic  Events 

Once  an  event  is  located  and  timed  by  an  automatic 
association/location  process,  the  next  step  is  to  access  time 
windows  containing  signals  from  the  event.  Data  within  signal 
time-windows  are  examined  to  determine  the  detection  status  of 
the  desired  signals.  Those  signal  windows  with  good  data  but 
without  detectable  signals  are  reduced  to  a  set  of  noise 
measurements.  Those  with  bad  data  are  flagged  as  inoperative. 
But  when  signals  are  detected,  the  system  extracts  the  signal 
waveform  for  present  and  future  quality  control  and  maintenance 
or  data  oases.  Also  detected  signals  are  reduced  to  a  set  of 
signal  measurements  for  the  on-going  real-time  operation  of 
the  VSC  system.  The  decision  tree  required  to  perform  this 
automatically  is  shown  in  Figure  1-3.  The  data  reduction 
accomplished  by  the  decision  tree  is  shown  in  the  left  margin 
of  the  figure. 


In  the  Identification  Experiment,  we  utilired  a  detector 
designed  by  Unger  (1973)  to  detect  and  time  short-period  signals. 
Unger  tested  his  detector  against  an  analyst's  detections 
finding  that  ninety  percent  of  the  signals  detected  by  both 
an  analyst  and  the  automatic  detector  were  timed  within  +_  0.5 
seconds  of  the  analyst's  pick.  Slowly  emergent  signals  were 
picked  several  seconds  late.  This  compared  favorably  for 
timing  with  the  more  conventional  Z  detector  of  Swindell  and 
Snell  (1977).  By  comparison,  the  Z  detector  times  signals 
1.70  seconds  late  on  an  average  with  a  standard  deviation  of 
-.Z3  seconds.  For  that  reason,  Urtger's  detector  was  selected 
tor  the  Event  Identification  Experiment  to  automatically  time 
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ar‘“  ecit  snort -period  ?  wave  signals  from  the  long  four 

»w5  .scores  initially  accessed.  Our  experience  indicated 
that  where  signals  were  detected,  they  were  accurately  timed 
*_..h  -nl\  negligible  false  alarms.  We  noted  however,  that 
impulsive  signals  much  shorter  than  the  selected  detection 
ga ^  a:id  slowiv  emergent  signals  were  often  missed  by  the 
aetector.  In  some  cases,  even  very  obvious  signals  were 
rnisaec  Decause  or  spikes  and  glitches  occurring  in  the  warm- 
-i0-se  Preceding  the  signal  window.  .Another  serious  problem 
was  also  observed  in  using  Unger’s  detector  for  automatically 
signals.  Extremely  grave  editing  errors  were  caused 
by  erroneously  selecting  mixed  events  and  multiple  phas es  . 

This  stemmed  from  the  lack  of  pos t- detection  validation  and 
a  too  simplistic  detection  strategy. 


•ve  consider  problems  caused  by  missed  signals  and  editing 
blunders  to  be  a  very  serious  flaw  in  our  present  software  for 
automatically  editing  signals.  We  also  believe  that  these 
problems  are  clearly  in  focus,  can  be  and  should  be  solved. 

That  goal  will  be  a  principal  objective  of  this  design  study. 

An  sxtsIls--ve  modirication  of  Unger's  detector  will  be  pre¬ 
sentee  in  a  later  section  of  this  report.  This  new  version 
of^ an  automatic  editor  should  correct  the  observed  faults  and 
De  recommended  for  future  implementation  into  the  VSC 
2>3wem.  One  important  lesson  learned  from  the  Event  Identifi¬ 
cation  Experiment  was  that  this  detector  was  an  important  factor 
-imiting  our  event  identification  performance. 
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SECTION  II 

^Eii_R.?TION  OF  UNGER'S  ANALYTIC  DETECTOR 


IHEORETICAL  3ACXGR0UND 


.  'he  baSiC  COnce?t  of  analytic  detector  is  to  represent 
.cloaca  as  the  real  part  of  a  sampled  complex  exponential 
--neuron.  3y  this  means  one  can  measure  signals  and  noise  as 
modulated  amplitude  and  phase  time  series.  For  example,  com¬ 
plex  seismic  sensor  noise  can  be  represented  as  a  sequence  of 
--  10m  ?hase  and  amplitude  modulations,  as  follows. 


nj  3  Nj  Exp  i<pj 


0  3  1,2,  ..  .  ,J)  . 


CII-1) 


Similarly  a  seismic  signal  occurs  within  at  1 
0"  points  on  the  record  as 


east  one  subset 


s.-  =  S.  Exp  i'ii. 

J  ]  3 

0 


(K< j  <  K+D) ,  and 
elsewhere 


"here  D  is  the  nominal 
measured  seismic  data, 


x . 


duration  of  the  seismic  signal.  The 
xj ,  is  tde  sum  of  signal-plus-noise. 

Sj  Exp  •  (IX- 33 


-he  detection  problem  is 
x,  and  .o  extract  and  measure 


to  time  the 
the  signal. 


start  of  the  signal, 
This  representation 


c  N  SCO 


:nc, 


.  1 


'U 


a'  a  s 


-ppi_sd  lO  the  detect  ton  and  measurement  of  seismic  sig- 
r.a.s  by  Farnbach  (1975)  and  Unger  [19  73).  One  of  the  earliest 
.reatments  or  this  method  was  given  by  Dugund j  i  [1958). 


-2 


*ne  application  of  this  type  of  phasor  representation  of 
gna_s  and  noise  is  practically  accomplished  by  applying  a 
5,_.ai  nicer  to  observed  data.  The  filter  passes  signal  and 
ise  energy  without  any  change  in  the  amplitude  but  with  a 
change  m  tne  phase  angle  spectrum.  This  operator,  called 
r.iibert  transform  is  described  by  Bracewell  [1565). 


ihis  operation  can  be  performed  efficiently  by  means  of 
a  fast  Fourier  Transform  (FFT) .  It  essentially  interchanges 
the  real  and  imaginary  parts  of  the  FFT.  This  is  followed  by 
inverting  the  data  spectrum  back  to  the  time  domain.  Alter¬ 
natively,  Quadrature  filtering  is  a  much  more  computationally 
erricient  procedure.  With  unit  amplitude  response  it  shifts 
each  rrequency  90°.  This  more  practical  operator  yields  re¬ 
sults  almost  identical  to  the  Hilbert  transform. 


Tne  time  series,  y^ ,  produced  by  the  Hilbert  transform 
sr  quadrature  filter  is  treated  as  the  imaginary  part  of  the 
Tne  data,  Xj ,  is  itself  considered  to  be  the  real  part 
cc  tne  complex  phasor  record,  I.,  where 


By  cnis  means,  the  record  of  data  is  transformed  into  an 
t rtnogonal  sequence  of  rectangular  complex  coordinates.  These 
-*o  uime  series  x.  and  y^  are  then  transformed  into  a  polar 
coordinate  representation. 
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The  tine  series  A.  is  a  measure  a£  the 
c£  the  data;  the  time  series,  ck  ,  the 
The  time  derivative  of  o.  is  a  measure 
dominant  frequency  of  seismic  data. 


instantaneous  envelope 
instantaneous  phase, 
of  the  instantaneous 


In  using  this  method  to  compute  the  instantaneous  fre¬ 
quency  ,  a  serious  stability  problem  was  pointed  out  by  Fambach 
and  Unger;  especially  when  the  envelope  is  small.  This  can  be 
seen  by  deriving  an  expression  for  the  instantaneous  frequency 
as  follows. 

r  _  1 


From  the  above  expression,  it  is  seen  that  estimates  of  f. 

3 

become  unstable  when  A.  becomes  very  small.  It  can  be  further 
seen  that  relative  variation  of  £  are  three  times  those  of  the 
envelope . 

5A 
A 

Thus,  I0?i  point-to-point  independent  fluctuations  of  the 
envelope  could  result  in  a  30%  error  in  frequency  estimates. 
This  problem  of  precise  frequency  measurement  must  be  solved 
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to  effectively  utilize  phase  measurements  in  the  design  of 
optimum  analytic  detectors.  This  will  be  discussed  in  a 
later  section. 


4 


3.  UNGER'S  CRITERIA  FOR  TIMING  SIGNALS 

Unger  f 1 9 7 S )  constructed  a  remarkably  simple  model  upon 
w.nich  to  base  his  analytic  detector.  He  analyzed  the  problem 
of  detecting  weak  signals  as  one  of  detecting  effects  produced 
by  the  interference  between  signal  and  noise  p’nasors  .  The  sig¬ 
nal  in  the  time  window  of  duration,  D,  is  represented  by  c 
fixed  signal  amplitude  level  and  initial  phase.  The  noise  is 
also  represented  as  a  fixed  level  but  with  random  point-by- 
point  transitions  of  phase. 

x j  =  N  Ex?  icfi .  +  S  Exp  iwAQ  T.  (K< j  <K+D)  CII-7) 

=  N  Exp  itjK  ,  elsewhere  on  the  record  (j  =1 , 2  ,  .  .  .  ,  J) 

The  fixed  phase  of  the  signal  is  arbitrarily  set  to  zero.  For 
the  purpose  of  analyzing  detections  this  is  without  any  loss 
of  generality.  Although  it  simply  models  the  problem  of  detec¬ 
ting  a  weak  signal  added  to  noise,  it  is  also  an  extreme  over¬ 
simplification  in  that  the  random  fluctuations  of  the  envelope 
are  neglected.  Also,  the  phase  of  short-period  seismic  signals 
is  generally  more  complex  than  a  single  simple  deterministic 
phase  modulated  pulse. 
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Even  more  interestingly,  if  the  signal  phase  changes  are 
unirorm  as  indicated  in  equation  (II-7)  or  indeed  can  be  pre¬ 
dicted  by  any  deterministic  model,  then  the  signal  can  be 
detected  as  a  stationary  initial  phase  condition.  Unger  gives 
the  probability  of  an  observed  stationary  phase  condition  in 
terms  of  measured  phase  fluctuations,  A<j>,  occurring  within 
prescribed  limits. 

1  -  i  COS-1  c|)  C0<|<1) 

?(-T  <  iOj  <T  )  * 

1  |<1  •  (II-9) 

The  expected  probability  of  meeting  the  above  detection 
conditions  under  noise  is  0.5  for  both  equations  (II- 8)  and 
(11-9} .  Under  these  equal  likelihood  of  detection  conditions, 
cetection  gain  advantage  of  6  dB  is  obtained  by  such  a  phase 
cetectcr,  provided  that  some  deterministic  phase  model  can  be 
obtained  for  short-period  seismic  signals .  Unger  (1973}  at¬ 
tempted  this  with  a  quadratic  phase  versus  time  model;  but 


1 


i 


found  that  seismic  signals  commonly 
fluctuations  not  too  dissimilar  from 
suit,  Unger  applied  equation  CH_8) 
ris  analytic  detector  which,  will  be 


demonstrated  random  phase 
seismic  noise.  As  a  re- 
as  the  model  for  design, 
described  in  the  next  sub 


section . 


C.  THE  DESIGN  OF  UNGER'S  ANALYTIC  DETECTOR 

As  a  baseline  for  utilicing  an  analytic  detector  in  an 
automated  edit  process,  we  will  describe  the  -design  of .Unger's 
detector  and  its  application  as  an  automatic  signal  editor  m 
the  Event  Identification  Experiment.  There,  the  detector 
package  was  inserted  in  a  supervisory  routine  CSPEED)  which 
directs  automatic  control  of  signal  timing  and  extraction  to  a 
detector  and  subsequently  to  a  set  of  signal  measurement 

algorithms . 

First  the  short -period  record  is  filtered  with  a  broadband 
filter  to  subdue  high  frequency  effects  such  as  spikes  and 
aliasing  and  low  frequency  effects  such  as  data  offsets  and 
drift.  This  was  done  without  distorting  signals  by  placing  the 
filter  cutoff  points  at  0 . 3  He  and  9.9  He  for  20  He  sampled 
data  and  at  0.3  He  and  4.9  He  for  10  He  sampled  data.  Jhis  ?re' 
filtering  was  followed  by  application  of  the  Unger  (19/8)  de¬ 
tector  . 

The  procedure  of  detecting  and  timing  the  onset  of  short - 
neriod  (S?)  signals  is  as  follows  (Figure  II-l).  First,  over 
a  specified  warm-up  period  (e.g.,  40  seconds),  the  peak  noise 
envelope ,  !n|  ,  is  established.  This  peak  envelope  is 

cosine  tapered  over  subsequent  waver:rm  points,  wun  a 
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'  ~“le  C3ns‘-ant  Ce  •  S  •  >  with  a  60 -second  time  constant, 
?ea.\  value  is  halved  at  50  seconds  and  equals 
-e.o  a  o0  seconds).  An  envelope  value  exceeding  the  tapered 
pea.s  vaiue  established  a  new  noise  peak,  unless  a  signal  de- 
v-c.ion  is  declared;  in  that  case  no  noise  peak  update  takes 
p.ace  until  the  signal  is  declared  to  be  terminated. 


.n  signa.  detection  is  called  whenever,  in  a  forward  look¬ 
ing  (.eadmg)  time  window  of  specified  length  (e.g.,  4  seconds), 
tne  probability  that  the  envelope  is  greater  than  the  tapered 


peak  noise  envelope,  P( j  rg Ct) | > | n | max)  exceeds  a  spedified 


-nresnoka,  TH1  (e.g.,  TH1  =  0.5).  When  this  probability 
reaches  its  maximum  the  algorithm  starts  looking  for  the  first 
signal  envelope  peak.  When  the  ratio  of  first  signal  envelope 
peax  and  tapered  noise  envelope  peak  exceeds  a  second  specified 
tnreshold,  the  SNR  threshold  TH2  (e.g.,  TH2  =  2  to  5  dB),  the 
signal  detection  is  confirmed  and  a  frequency-dependent  step- 
oack  is  performed  to  determine  the  signal  onset  time. 


The  stepback  procedure  (Figure  1 1 - 2)  is  based  on  the  ob- 
se./a.^on  that  in  most  cases  the  first  signal  envelope  peak 
.at  tj  occurs  within  one  signal  period,  and  frequently  at 


approximately  5/4  period,  after  the  signal  onset  (at  tn) .  I 


n 


a  high -SNR  waveform  the  signal  onset  time  is  most  accurately 
-wund  by  detecting  the  rirst  maximum  or  minimum  of  the  signal's 
ins cantaneous  value  (at  t.) ,  and  stepping  back  1/4  period 
(  =  0  .25/ instantaneous  frequency  at  t.) .  For  low- SNR  waveforms 
tne  rirst  quarter  period  may  be  obscured  by  noise;  in  that 
w-se  *e  Seep  back  5/4  mean  period  (=0 . 75/mean  frequency  at  t^) 
rrom  tne  rirst  signal  envelope  peak  at  td.  The  mean  frequency 
v..e  c-osed-.orm  derivative  of  the  phase  regression  polynomial 
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is  star,  ted  a:  t„ ,  i.s. 
first  quarter  period  i 
exceeds ,  by  a  third  th 
immediately  preceding 


The  search  for  the  first  quarter  period 
,  at  0.3  mean  period  before  t  ^ ;  the 
s  detected  when  its  maximum  or  minimum 
reshold,  TH5  C e.g.,  TH5  =  1  dB)  ,  the 
noise  in  the  one-second  time  interval 


)  • 


--  second  threshold  (the  SNR  threshold)  is  not 
ried,  the  detection  is  annulled  and  the  noise  peak  value 
updated  with  what  at  rirst  was  believed  to  be  the  signal 
envelope  peak.  Thereafter,  the  noise  peak  is  updated  as 
'until  the  next  supposed  signal  detection,  etc. 


satis  - 
is 

usual 


t 


The  signal  end  time  is  found  as  the  moment  of  the  first 

envelope  minimum  occurring  either  after  P(|r  (t)|>|n!  )  falls 

.  .  .  _  s'-'  max' 

oeiow  its  threshold,  or  after  the  signal  duration  exceeds  a 

specified  maximum,  whichever  is  first.  If  this  envelope  is 

updated  with  this  value,  the  noise  peak  updating  and  signal 

detection  resume  as  normal.  In  principle  this  procedure 

enables  the  detection  and  timing  of  later  phases  and  other 

signals  in  the  coda. 


r rom  the  previously  described  model,  several  important 
advantages  can  be  cited.  For  signals  less  than  6  dB  over 

peaks  oi  estimates  oi  the  P ( A ^ > N )  correspond  to  maxima 
ci  5/N<2.  Thus,  these  indicators  of”weah  signals  satisfy  the 
maximum  likelihood  condition  independently  of  the  statistical 
■‘Oise  distribution,  .-or  signals  greater  than  6  dB  over  the 
noise,  >-ne  detector  saturates  .  Thus  ,  the  detector  orovides 
a  robust  means  of  detecting  and  initially  timing  the  larger 
s^gnais.  In  principal,  this  detector  should  optimally  time 
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:he  arrival  of  short -period  P  waves.  It  should  be  robust, 
a.so ,  because  an  estimating  the  probability  of  the  envelope 
exceeding  noise,  it  counts  the  number  of  times  that  the  sig¬ 
nal  exceeds  noise  in  a  leading  time  window  containing  the 
signal.  If  properly  implemented,  it  is  potentially  insensitive 
to  unusually  large  and  occasionally  erratic  envelope  fluctua¬ 
tions  m  the  noise  or  signal  windows.  3y  careful  design,  it 
should  be  insensitive  to  spikes,  glitches,  and  other  large 
amp  1 i tude  errors . 


D.  POST-MORTEM  EVALUATION  OF  UNGER'S  ALGORITHM  AS  AN  AUTOMATIC 
*  SIGNAL  EDITOR  . '  ‘  " 


Our  analysis  of  errors  in  applying  Unger's  automatic  de¬ 
tector  as  a  short-period  signal  editor  lead  us  to  suggest  the 
tallowing  requirements  for  acceptable  performance. 

•  Maintain  a  specified  acceptable  false  alarm  rate  in 
retrieving  signals  associated  with  a  known  event. 

•  Detect,  time,  and  validate  almost  as  many  signals 
as  a  seismic  analyst. 

•  Separate  signals  of  a  known  event  from  mixed  signals 
of  other  events  and  other  seismic  phases. 


Of  these  minimum  requirements,  we  can  only  claim  that  the 
tirst  one  was  met  by  our  Event  Identification  Experiment  ex¬ 
perience.  Thererore,  the  above  requirements  for  automatic 
editing  necessitate  a  re-design  of  the  analytic  detector 
algorithm.  The  following  brief  post-mortem  evaluation  is 
pro  viced  as  background  for  the  re-design. 


Eger's  analytic  detector  produced  negligible  false 

;::7"S  a‘?d  a;n°St  a11  0f  tiie  validated  detected  signals  were 
:;;crate.y  times  including  many  cases  of  barely  visible  sig- 

‘‘rrl'  ‘^US’  the  first  re<lui*ement  of  an  acceptable  false 
^-arn  rate  and  part  of  the  second  requirement  of  accurate 
ttmtng  was  satisfied  by  Unger's  analytic  detector. 

Miss ea  signals  caused  a  very  serious  problem  in  identify- 
"*3  *ea,  events.  The  detector  algorithm  operated  onlv  on 

en\  eiope  measurements  of  the  seis-Hr  nor*  ,  t, 

.„,<C4.1v  seism.c  aata,  wnereas  ,  an  analyst 

:c::  :y  ,per*'“  °» *  -»<*  <>1™™.  set  oi  sigatl  char. 

:  etasttas.  Given  .  reasonable  false  ala™  rate  we  estate 

.  Unger  s  aetsctor  could  only  detect  at  a  level  of  705  of 
-  -meal  analyst  detection  capability.  For  this  reason, 

*  recommend  that  detection  of  signals  be  based  on  a  more 

‘■V!‘SVet  0t  mUltiVariate  “Plitude  and  frequency  measure- 
men^s  Oi  ground  motion. 

...  "e  °bae"ed  “”her  reason  why  some  signals  obviously 
S.t0.an  analyst  were  missed  by  the  automatic  detector. 

detec:0r  al*orith-  mimiced  as  an  unknown 

P  -tsnsient  added  to  stationary  Gaussian  noise.  On 

;.tat  oasts  an  envelope  threshold  is  set  to  detect  signals  and 

*  '  ,a*aras'  Glvea  such  »  ”od«l.  such  an  optimum 
“  "  '  ,  °  be  S*able  and  the  operating  characteristics 

u  a  then  be  determinable.  In  practice  often  such  idealised 

d,’.  ?  arS  “h  tTU#-  F"  «“*»•  very  large  amplitude 

"s  cccas  tonally  occur  from  environmental  effects  such 
*?  ”S  °r  f'°”  e-sc-ronic  malfunctions.  Because  Unger's 

:'s°  ”  gaUges  Zne  31axi-'"U"’  envelope  of  observed  noise  to 

-e-ermcne  tne  nominal  noise  level  such  large  glitches  or 


s?i.<es  can  serve  to  effectively  shut-off  the  detector.  In 
such  cases  only  very  large  events  can  be  detected.  Only  sig¬ 
nals  much  larger  than  the  noise  glitch  would  be  detected. 

"""S  nappenea  to  us  enough  times  in  the  Event  Identification 
experiment  to  recommend  modifying  Unger's  procedure  of  gauging 

the  interfering  noise  level  and  replacing  it  by  a  more  robust 
procecure . 

-ven  when  the  noise  data  were  of  excellent  quality  we 
encountered  numerous  cases  of  missed  visit  '  signals.  This 
was  one  of  many  problems  caused  by  the  diverse  nature  of  seis- 
•tiic  signals,  m  this  case  by  emergent  seismic  signals.  Unger's 
cetector  is  designed  to  be  optimum  when  signal  starts  at  some 
point  K^and  maintains  a  fixed  level  for  duration  D  and  then 
snuts  orf  or  gradually  decays.  This  assumption  is  not  at  all 
optimum  tor  detecting  emergent  seismic  signals.  In  that  case, 
Unger's  maximum  noise  level  estimate  will  be  sequentially  up¬ 
dated  when  the  emergent  signal  rises  above  the  noise  but  re¬ 
mains^  below  the  detection  threshold  of  a  signal.  As  a  result, 
the  threshold  is  continuously  raised.  A  gradually  emerging 
signal  is  then  either  missed  altogether  or  detected  with  a 
•  e.,.  _arge  time  delay.  This  is  another  example  where  Unger's 
noise  level  updating  procedure  is  not  sufficiently  robust  to 
-e.ee.  emergent  signals.  It  is  recommended  that  the  presently 
appi.ed  noise  level  updating  procedure  be  replaced  by  a  more 
-coust  noise  estimation  procedure  which  is  less  sensitive  to 
gl Irenes  and  spikes  and  to  interference  effects  produced  at 
t.te  beginning  of  emergent  signals. 

-n^this  regard,  we  point  out  that  Unger's  algorithm  was 
w..g^na_ly  designed  as  a  tront-end  continuous  waveform 
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SECTION  III 

ANALYSIS  OF  AN  ADVANCED  ANALYTIC  DETECTOR 


In  the  preceding  section,  our  post-mortem  evaluation  of 
Unger's  cetector  indicated  excellent  performance  in  tim:ng  sig¬ 
nals  .  That  is  for  signals  which  both  the  automatic  detector 
and  analyst  detected.  However,  in  the  case  of  emergent  com¬ 
plex  seismic  events,  accurate  timing  presented  a  problem  which 
needs  to  be  solved.  But  the  real  problem  with  the  detector  is 
the  missed  signal  problem. 

The  automatic  detector  misses  some  large  signals  obvious 
to  an  analyst  and  many  more  small  signals  which  can  be  detec¬ 
ted  by  a  seismic  analyst.  Some  of  the  more  obvious  automatic 
detector  misses  could  have  been  avoided  by  minor  modifications 
to  the  detector. 

The  severity  of  the  problem  was  such  that  for  about  40% 
of  the  events,  no  signals  were  detected  by  the  automatic  de¬ 
tector.  In  those  cases,  discriminants  had  to  be  derived  by 
networking  noise  estimates.  This  was  to  put  upper  limits  on 
magnitude  measurements  which  were  derived  to  compute  discrim¬ 
inants.  Although  this  approach  produced  good  results,  it  be¬ 
came  obvious  that  signals  missed  by  the  automatic  detector 
were  primal  m  limiting  event  identification  performance. 
Clearly,  more  than  minor  modifications  to  Unger's  detectoi 
would  be  required  to  achieve  a  satisfactory  level  of  perfor¬ 
mance  . 
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Of  the  signals  examined  in  the  Event  Identification 
Experiment,  it  is  estimated  that  10 %  of  them  which  were 
missed  by  the  automatic  det<  ctor  were  obviously  large  sig¬ 
nals  easily  detected  by  an  ;<ialyst.  These  signals  probably 
could  have  been  detected  automatically  by  relatively  minor 
mcdif icaticns  of  Unger’s  detector.  Suggested  modifications 
are  as  follows : 

•  Use  ordered  statistics  to  gauge  noise  thereby 
avoiding  influence  of  occasional  large  glitches 
and  spikes. 

•  Combine  Unger's  detector  with  a  conventional  I 
detector  to  detect  impulsive  earthquakes  and  ex¬ 
plosions  . 

•  To  avoid  missing  or  late  timing  of  emergent  signals, 
count  strings  of  ascending  envelope  measurements  as 
above  threshold  if  terminated  by  at  least  one  en¬ 
velope  value  above  the  threshold. 


Even  with  these  changes  in  Unger's  detector  algorithm  we 
would  still  expect  to  miss  about  20%  to  50%  of  the  events 
which  could  be  detected  by  an  experienced  seismic  analyst. 
This  follows  from  the  fact  that  Unger's  algorithm  detects  on 
sustained  signal  power  over  some  nominal  signal  duration,  D. 
.is  such  it  is  no  different  than  most  seismic  detectors  which 
attempt  to  detect  signals  solely  by  their  enhanced  broadband 
power.  At  reasonable  false  alarm  rates,  such  detectors  can¬ 
not  be  expected  to  detect  more  than  70%  to  30%  of  the  low 
signai-to-noise  rate  signals  detected  by  seismic  analysts, 
Swindell  and  Snell  (1977] . 
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bridge  this  gap  of  automatically  detecting  weak  sig- 
h  the  something  approaching  the  efficiency  of  a  seis- 
yst,  we  will  analyte  the  design  of  a  new  multivariate 
detector.  This  will  be  discussed  in  the  following 
ons  . 


A.  GENERALIZED  ANALYTIC  SIGNAL  AND  NOISE  MODEL 

The  theoretical  background  for  the  design  of  analytic 
detectors,  described  in  Section  II,  is  generalized  to  encom¬ 
pass  tne  diverse  types  of  seismic  signals  expected  from  earth¬ 
quakes  and  explosions.  These  include  multiple  transmissions 
cue  to  a  complex  source  or  to  different  propagation  paths, 
variable  length  duration  of  signals,  multiple  pulse  content 
or  signals,  and  the  fading  coda  scattering  associated  with 
signals . 

Ambient  noise  preceding  the  time  window  searched  for  sig¬ 
nals  is  generalized  to  possibly  include  spikes  and  glitches 
due  to  intermittent  system  malfunctions ,  environmental  effects 
at  the  receiver  site,  and  local  seismicity. 

1.  .Ambient  Noise  Model 


each  record  accessed  by  the  automatic  signal  editor  will 
contain  at  least  two  minutes  of  short-period  data  to  be 
scorched  ror  signals.  One  half  to  one  minute  of  noise  data 
will  re  provided  to  gauge  the  noise  statistical  distribution. 
..*ese  noise  statistics  will  be  presumed  to  be  fixed  over  the 
balance  of  the  record  possibly  containing  the  desired  signal, 
equation  C *  - ~ 1 3  describes  the  complex  seismic  sensor  noise  at 
-.ie  3  tat.on.  »his  ambient  noise  model  is  generalized  to 
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include  spikes,  glitches,  and  local  events.  These  will  oc¬ 
cur  as  short  bursts  of  energy  which  complicate  the  signal 
detection  decision  function.  The  noise  model  is  generalized 
to  include  large  non-stationary  energy  fluctuations  in  the 
1  noise.  Per  this  we  obtain 

'  kN 

n,  »  CN-  +  Z  Go)  Exp  io,  ( i  =  1 , 2  ,  .  .  .  ,  J)  (III-l) 
-  £»k.  1  3  L  J 

l  A 

,  where 

<5  *1  for  m=0  or  zero  elsewhere;  and 
m 

<  k., ,  k, , . . . , are  N  random  energy  spike  occurrences. 

Unger  73)  performed  a  study  of  the  distribution  of  Nj 
for  seismic  noise.  Normalizing  the  noise  fluctuations  by 
j  dividing  by  the  standard  deviation,  he  interpreted  the  seismic 

i  noise  trace  as  a  stationary  Gaussian  distribution.  In  that 

case,  the  envelope  of  that  noise  would  be  a  Rayleigh  distribu¬ 
tion.  Since  measurements  of  signal  envelope  are  primarily 
used  to  derive  magnitudes,  he  also  determined  the  distribu¬ 
tion  w  =  log10>L;  the  noise  magnitude  distribution. 

Based  on  the  analysis  of  one  hour  of  seismic  noise  at  the 
Korean  Seismic  Research  Station  (XSRS)  broken  down  into  102.4 
second  data  segments,  Figure  1 1 1  - 1  shows  Unger's  results  in 
fitting  measured  noise  magnitude  distributions  to  the  distri¬ 
bution  derived  from  a  stationary  Gaussian  noise  assumption. 
Unger  tested  the  significance  of  the  fit  of  measured  noise  to 
a  Rayleigh  distribution.  At  the  53  significance  level  303  of 
the  measured  distributions  were  consistent  with  the  Gaussian 


HNS CO, 


:xc. 


Ill -4 


noise  hypothesis.  Of  the  20$  of  rejections,  the  Gaussian 
hypothesis  slightly  underestimated  the  number  of  large  magni¬ 
tude  fluctuations.  Thus,  these  rejected  samples  appeared  as 
almost  Gaussian  distributed. 

Based  on  the  Gaussian  noise  assumption,  certain  facts 
can  be  ascertained  about  magnitude  measurements  of  seismic 
noise.  Some  of  these  were  given  by  Unger  (1978)  others  were 
derived  from  graphical  analysis  of  his  results.  The  most 
probable  occurrence  of  magnitude  is  at  the  log.,  q  (RMS)  of  seis¬ 
mic  noise.  The  standard  deviation  of  large  noise  fluctuations 
above  the  log. q (RMS)  of  noise  is  approximately  0.20  magnitude 
units.  Measurements  at  seven  Seismic  Research  Observatory 
(SRO)  stations  (Strauss  and  Weltman,  1977;  Weltman,  et  al., 
1979)  indicated  that  peak  one-second  noise  amplitude  measure¬ 
ments  had  estimated  magnitude  standard  deviations  between 
0.17  and  0.23;  and  four  stations,  between  0.13  and  0.16. 

The  shape  of  the  distribution  of  seismic  magnitudes 
based  on  Unger's  observations  as  well  as  on  the  Gaussian 
noise  assumption  is  highly  skewed.  Only  one  third  of  the 
envelope  measurements  are  expected  above  the  RMS  noise  (the 
most  probable  noise  magnitude) .  Of  the  two  thirds  of  the 
envelope  measurements  below  the  RMS  noise  level,  the  RMS  of 
these  negative  magnitude  fluctuations  is  approximately  0.5 
magnitude  units.  On  this  basis  it  is  estimated  that  in  20$ 
or  a  time  window  containing  the  seismic  noise  envelope  mea¬ 
surements  are  lower  by  more  than  a  half  magnitude  than  the 
r.cise  RMS  level. 

At  .intervals  of  time  where  such  low  level  noise  envelopes 
are  encountered,  a  skilled  analyst  can  possibly  recognise  and 
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5-5..J..5  by  -n uermittent  changes  in  the  character,  such 
as  by  changes  of  frequency  and  the  duration  of  such  changes. 
Perhaps  this  night  explain,  in  part,  the  20%  of  signals  which 
v--..  oe  ..erected  by  skilled  analysts;  but  which  are  always 
mss ed  by  detector  algorithms  based  solely  on  the  signal 

-w.wl.  ii  suggests  that  weak  signals  might  sometimes 
be  detected  at  noisy  stations  at  levels  as  much  as  one  half 
magnitude  Higher  than  the  expected  signal  magnitude.  If  fur- 
'"“s*  suggests  that  the  algorithm  for  automatically  editing  and 
measuring  signals  be  based  on  more  than  energy  level  criteria 
alone . 


As  for  our  automated  signal  editing  procedure,  the 
generalication  of  the  noise  model  of  equation  Q I I I - 1 )  can  be 
applied  to  the  following  tasks. 

•  Design  robust  noise  parameter  estimation  techniques 
which  are  insensitive  to  occurrences  of  spikes, 
glitches,  and  local  events. 

•  Automatically  detect  and  flag  malfunctioning  sta¬ 
tions  . 


The  cirst  step  in  this  analysis  is  to  access  the  one 
minute  window  containing  seismic  noise.  This  noise  data 
precedes  the  t.nree  minute  window  possibly  containing  signal 
anc  coda.  As  is  usual  in  robust  estimation  procedures, 
noise  measurements  are  analyted  by  application  of  ordered 
a .a-.a wics .  This  is  to  avoid  influencing  the  result  by 
neavily  weighting  large  deviations  due  to  spikes,  glitches, 
arc  local  seismic  events,  as  shown  on  equation  (III-l). 
envelope  measurements  of  the  noise  are  sorted  from  smallest 
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tc  -argest.  By  raxing  the  logarithm  to  base  ten  of  the  noise 
envelope  measurements,  the  envelope  measurements  are  trans - 
rorned  into  a  set  of  ordered  noise  magnitude  measurements. 
Applying  Unger's  derived  distribution  for  log-envelopes, 
the  following  percentiles  are  determined  by  interpolating  the 
crcered  set  of  measured  noise  magnitudes. 

•  ihe  66  percentile  level  of  smaller  magnitudes  cor¬ 
responds  to  log]_g(RMS)  of  the  noise  trace  and  cor¬ 
responds  to  the  most  probable  occurring  envelope 
magnitude . 

•  ihe  S3  percentile  level  minus  the  66  percentile 
corresponds  to  the  standard  deviation  of  positive 
noise  magnitude  fluctuations,  expected  to  be  ap¬ 
proximately  0.20  magnitude  units. 

•  The  66  percentile  minus  the  20  percentile  level 
corresponds  to  the  standard  deviation  of  negative 
noise  magnitude  fluctuations,  expected  to  be  ap¬ 
proximately  0.50  magnitude  units  for  Gaussian 
noise . 

The  expected  values  of  the  positive  and  negative  fluctua¬ 
tions  of  noise  magnitude  were  observed  by  Unger  at  KSRS,  and 
correspond  to  what  is  expected  for  ideal  bandpass  filtered 
Gaussian  noise.  3y  contrast,  noise  from  ocean  generated 
microseisms  is  characterized  by  a  single  sharp  spectral  peak 
between  0.15  and  0.55  He.  Rice  (1954)  analyzed  this  situation 
zr  a  periodic  process  added  to  random  Gaussian  noise.  If  the 
periodic  process  is  large,  the  noise  distribution  of  envelope 
measurements  shifts  from  the  Rayleigh  distribution  (the  envel- 
:pe  distribution  of  Gaussian  noise)  to  a  normal  distribution 
,v - -h  standard  deviation  of  the  RMS  of  the  additive  Gaussian 
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noise.  The  mean  of  the  normal  distribution  is  centered  at 
me  amplitude  of  the  periodic  process.  This  is  shown  in 
ngure  Thus,  unless  the  influence  of  the  micros eismic 

pea.<  is  minimizes  by  filtering  or  some  other  means,  Unger's 
app .ication  of  the  Guassian  noise  model  would  have  to  be 
modelled,  in  the  case  of  a  large  nearly  periodic  component 
m  the  seismic  noise,  noise  magnitude  fluctuations  would  be 
nearly  symmetrical  ror  positive  and  negative  deviations  and 
a. most  normally  distributed.  For  a  spectral  peak  at  three 
times  the  RMS  Gaussian  noise,  the  standard  deviation  of 
envelope  magnitudes  would  be  about  0.15  magnitude  units.  The 
oRC  noise  magnitude  measurements  suggest  that  at  least  some 
seismic  stations  would  be  close  to  this  alternative  noise 
model.  In  those  cases,  techniques  would  need  to  be  applied 
to  minimize  the  effect  of  sharp  spectral  peaks  occurring  in 
the  noise  to  optimize  detection  performance. 

A  systematic  procedure  for  applying  ordered  statistics 
to  magnitude  measurements  is  given  by  Sax,  et  al.  (1979)  and 
is  illustrated  in  Figure  1 1 1  -  3  .  Based  on  Unger's  envelope 
study,  the  percentiles  applied  to  noise  data  which  is  nearlv 
oaussian  are  given  above.  By  applying  this  procedure,  magni¬ 
tude  measurements  of  seismic  signals  and  noise  are  reduced 
approximately  to  homogeneous  normal  statistics.  The  dis¬ 
crimination  operating  results,  after  using  this  normaliza- 
wicn  procedure  shown  on  Figure  1 1 1  -  3 ,  could  be  closely  pre- 
cictec  oy  application  or  normal  error  theory.  For  this  rea- 
;wn,  it  should  be  applied  to  the  problem  of  detecting  sig¬ 
nals,  where  the  distribution  of  noise  and  signal  magnitudes 
may  be  highly  skewed. 
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“‘1S  ordered  statistical  analysis  can  be  used  to  obtain 
..te  nest  probable  occurrence  of  noise  measurements  and  the 
'  tandard  deviation  of  positive  or  negative  fluctuations . 

This  reduces  deviations  from  the  most  probable  noise  oc¬ 
currence^  unit  normal  statistics.  In  singal  windows  a  set 
or  s igniricant  deviations  from  a  unit  normal  distribution 
can  be  used  to  detect  and  time  signals. 

Unusual  deviations  can  also  detect  and  diagnose  auto- 
ma._v.aily  those  stations  which  are  malfunctioning.  Note  that 
unusually  large  envelope  fluctuations,  Gk,  in  equation  CIII-1) 
can  oe  seen  as  unusually  large  noise  magnitude  deviations, 

but  will  not  influence  the  noise  parameters  used  to  detect 
signals . 

in  validating  a  properly  functioning  station,  the  median, 
positive,  and  negative  standard  deviations  can  be  checked 
against  historical  data.  If  these  parameters  fall  within 
acceptable  limits,  e.g.,  falls  within  a  99%  acceptance  re¬ 
gion,  the  station  noise  data  are  further  checked  for  large 
intermittent  errors,  .  If  the  average  of  observations 
above  the  one  standard  deviation  limits  is  more  than  two 
standard  deviations,  then  those  measurements  closest  to  the 
outner  average  are  counted  as  noise  spikes  or  glitches,  G.  . 
Positive  Gk's  may  indicate  spikes  or  clipped  data.  Negative 
Gk ' s  indicate  :ero  returns  or  intermittent  dead  spots. 

-r  tne  numDer  of  G^s  exceed  an  acceptable  value,  the  station 
is  i__ass  ined  as  a  malfunction. 

As  part  of  the  quality  control  built  into  the  automatic 
'-mg  detector,  the  noise  window  is  searched  for  soarselv 
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occurring  sp 
spikes  occur 
terpreted  as 


i.<es,  G,_..  I:  more  than  a  specified  number  of 
above  three  standard  deviations ,  these  are  in- 
sparse  occurrences  of  spikes.  In  that  case, 


one  amplitude  and  phase  measurements  of  the  record  are 


smoothed  (similar  to  ' liftering'  in  cepstral  analysis]  to 
minimi ce  the  effect  of  sparse  occurrences  of  .  This  will 
usually  have  small  distortional  effect  on  signals  since 
-heir  duration  will  be  at  least  one  half  second  or  greater. 

-he  same  time  it  will  minimice  the  influence  of  sparsely 
occurring  spikes  or  much  smaller  duration. 


The  above  quality  control  criteria  is  merely  a  prelim¬ 
inary  specirication  of  the  procedures  needed  to  validate 
properly  functioning  stations  and  to  condition  data  collected 
by  marginally  runctioning  stations.  Obviously,  these  criteria 
neec  to  be  tested  on  stations  exhibiting  various  types  of 
malfunctions .  Such  tests  will  result  in  modified  optimized 
procedures  for  validating  stations  and  conditioning  accessed 
signal  records.  In  the  Event  Identification  Experiment  this 
quality  control  was  done  by  human  analysts.  In  the  VSC  sys¬ 
tem,  we  anticipate  that  the  station  editing  quality  control 
will  need  to  be  done  automatically  in  order  to  quickly  re- 
cuce  the  volume  or  seismic  data  to  manageable  proportion. 

Generalized  Signal  Model 

Having  passed  through  noise  quality  control,  signal  data 
are  processed  to  detect  and  time  seismic  signals.  In  order 
to  obtain  valid  signal  measurements,  such  detections  need  to 
be  modeled  with  sufficient  generality  to  verify  that  mea¬ 
surements  will  be  made  on  proper  P  wave  signals.  The 
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generalized  signal  no del  must  therefore  include  effects 
produced  by  multiple  phase  reception  of  the  signal,  scattered 
coda  characteristics ,  and  instrument  malfunctions. 

The  multiple  signals  accounted  for  by  the  model  should 
include  the  following. 

•  Complex  source  signals  and  depth-phases,  each 
propagating  world-wide  as  first-motion  P  waves 

•  Surface-  and  core-reflected  phases  propagating 
world-wide  as  later  phases  of  each  P  wave  signal 

•  Long-distance  propagated  PKP  and  PKIXP  core  phases 

•  Localized  multiple  transmissions  and  receiver- 
scattered  phases  not  correlated  at  other  stations 
of  the  network 

•  Mixed  signals  from  another  event. 

This  model  describing  signal  phases  which  might  be  detected 
in  signal  time  windows  represents  a  complicated  situation. 
Therefore,  the  signal  model  is  changed  to  handle  this  situa¬ 
tion  . 

The  simpler  single  phase  signal  model  of  equation  (II-2) 
is  generalized  to  that  of  multiple  phase  transmissions  from 
a  complex  source.  As  in  the  noise  model,  we  also  include 
to  represent  intermittent  system  noise  in  the  accessed  signal 
window.  The  generalized  signal  model  is  given  by  equation 
(HI-2)  . 
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.he  model  is  generalized  to  MN  signal  source-transmissions 
propagated  to  each  seismic  station.  This  model  indicates  a 
complex  source  generating  a  sequence  M  signals  over  a  time 
interval  of  one  or  two  minutes.  Each  signal  is  propagated 
as  a  primary  P  wave  (n=  1)  followed  by  a  sequence  of  later 
arriving  secondary  phases  (n=2,...,N),  provided  that  phases 
such  as  ?P,  PcP,  etc.  occur  within  the  time  frame  of  the  ac¬ 
cessed  record.  Thus,  the  sequence  of  start-times,  K  ,  can 
•  •  mn 

be  viewed  in  a  network  sense  as  a  sequence  of  origin  times, 
each  augmented  by  the  propagation  time  delay  of  P  waves  or 
secondary  phases. 


Complex  explosion  sources  are  realized  by  detonating  a 
sequence  of  explosions  at  the  source.  Complex  earthquakes 
are  commonly  observed  from  some  seismic  regions.  Sax  (1979) 
ooserved  such  sequences,  which  from  array  processing  were 
apparently  propagating  as  source  delayed  P  waves  in  accord¬ 
ance  with  the  above  model.  The  duration  of  signals,  D,  is 
also  viewed  as  a  source- dependent  parameter  which  varies 
nominally  between  one  half  to  ten  seconds.  Small,  high 
stress  and  high  stress  drop  sources  woulc.  tend  toward  the 
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lew  end  of  the  duration  range  as  simple  and  highly  impul¬ 
sive  events.  Complex  signal  waveforms,  associated  with  large 
low  stress  drop  events  from  highly  heterogeneous  source  re¬ 
gions  ,  would  tend  toward  the  high  end  of  the  duration  range. 
The  design  of  the  detector  must  be  sufficiently  general  to 
cope  with  these  realistic  but  more  complex  signal  editing 
situations.  The  situation  is  further  complicated  in  that 
signal  magnitudes  associated  with  a  complex  source  may  also 
vary  considerably,  i.e.,  the  magnitude  of  later  source  emis¬ 
sions  may  be  larger  than  earlier  ones. 

This  obviously  complicates  the  detection  validation 
problem  of  correctly  associating  independent  edits  at  one 
station  to  those  at  other  stations  of  the  network.  This 
problem  of  validating  station  detections  is  considerably 
simplified  by  using  the  criteria  of  a  maximum  acceptable 
travel  time  anomaly  between  the  arrival  of  an  elemental 
source  phase  at  any  two  stations  of  the  network.  It  is  im¬ 
portant  to  validate  signal  edits  not  only  to  assure  propaga¬ 
tion  characteristics  consistent  with  the  given  event  location 
but  also  provides  a  sound  relatively  easy  basis  for  auto¬ 
matically  shifting  the  event  focus  to  a  position  which  mini¬ 
mises  observed  travel  time  anomalies. 

In  addition  to  utilising  this  criteria  of  propagation 
consistency  to  validate  automatically  edited  signals,  we  will 
also  utilise  magnitude  consistency.  The  difference  in  magni¬ 
tude  of  network  associated  phases  should  be  maintained  within 
a  maximum  acceptable  tolerance.  Although  source  magnitude 
consistency  is  a  weaker  criteria  than  propagation  consistency 
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che  variability  of  magnitudes  between  stations  , 
important  to  utilize  it  to  avoid  large  blunders 
raise  associations  of  an  event  with  large  receiver- 
phases  or  mixed  signals. 


To  complete  the  analytical  signal  model  we  include  a 
term  to  encompass  other  apparent  signals  which  can  be  de¬ 
tected  at  a  station.  These  signals  do  not  correlate  across 
the  network  and  therefore  cannot  be  validated  by  means  of 
propagation  consistency  criteria.  In  some  cases,  signal 
measurements  of  these  uncorrelated  phases  could  be  very 
misleading . 


These  uncorrelated  signals  from  the  source  can  in  some 
cases  be  associated  with  the  desired  event.  There  could  be 
one  or  two  delayed  secondary  images  of  the  source  caused  by 
rapid  upper  mantle  increases  in  the  propagation  velocity. 
Such  phases  can  be  large  but  would  be  seen  only  in  a  narrow 
distance  band  from  the  source  and  therefore  would  not  be 
generally  correlated  across  the  network.  Also,  strong  sig¬ 
nals  could  be  produced  from  an  inhomogeneous  source  medium 
as  higher  order  multipole  components.  These  would  only  be 
seen  m  narrow  azimuth  bands  and  also  not  generally  be  cor¬ 
related  across  the  network. 


Other  uncorrelated  signals  are  not  directly  related  to 
••‘2  source  and  can  be  caused  by  large  envelope  fluctuations 
cue  to  receiver  scattering  by  interfering  signals  from  some 
otner  known  event,  singular  observations  of  an  unknown  event, 
i.e.,  local  or  regional,  or  noise  false  alarms.  In  the 
latter  case,  we  may  have  two  events  with  overlapping  records 
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at  stations.  It  is  important  that  these  signals  from 

..L.xe--  events  ce  properly  sorted  cut  by  0Ur  post  detection 
a  n  a  1  v  s  1 5  , 


As  for  representing  these  detectable  signals  which  oc¬ 
cur  singularly  and  cannot  be  verified  as  propagation  con¬ 
sistent,  we  will  generalize  the  signal  model  of  equation 
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he  start  times  Kn  are  of  false  signals.  These  cannot  be 
venned  either  as  propagation  consistent  or  as  source  mag¬ 
nitude  consistent.  They  will  generally  be  rejected  as  noise 
One  or  the  advantages  of  this  model  is  that  the  decision 
threshold  can  be  set  much  lower  to  detect  weak  signals.  Bv 
^Plying  post  detection  analysis,  false  alarms  from  ambient 
noise  and  other  undesireable  detections  can  be  weeded  out  of 
-ne  automated  signal  editing  process.  The  purpose  of  rep- 
--sen-.ng  raise  signals  of  type  Uv  in  equation  CHI-3),  is 
to  establish  criteria  for  identifying  such  signals  bv  the 
lacs  of  network  propagation  consistency,  source  magnitude 
consistency,  or  other  criteria  and  to  keep  such  false  detec¬ 
tions  rrom  contaminating  the  signal  measurements  to  be  used 
ror  event  discrimination. 
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.AMPLITUDE  AND  PHASE  STATISTICS 
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Multivariate  statistics  are  needed  to  effectively  re¬ 
trieve  small  signals  below  the  RMS  level  of  seismic  noise. 
This  need  stems  from  the  often  observed  gap  between  the  sig¬ 
nal  retrieval  capability  of  a  seismic  analyst  compared  to 
that  of  an  automatic  algorithm  which  is  based  solely  on  power 
fluctuations . 


In  using  Unger's  algorithm  to  time,  measure,  and  retrieve 
signals  in  the  Event  Identification  Experiment,  we  failed  to 
retrieve  one  or  more  signals  for  about  40%  of  the  events 
analyzed.  In  several  percent  of  the  cases,  we  not  only  missed 
retrieving  signals  of  a  desired  event  but  committed  the  more 
serious  error  of  retrieving  much  larger  signals  from  some 
other  unknown  event  occurring  on  the  same  seismic  record. 

About  25%  of  these  missed  signals  and  most  of  the  serious 
editing  blunders  could  have  been  avoided  by  modifying  Unger's 
algorithm;  in  effect  using  more  robust  techniques  to  imple¬ 
ment  the  algorithm  and  performing  post  -  detection  analysis  to 
validate  signals.  Yet,  at  reasonable  false  alarm  rates,  ap¬ 
proximately  50%  of  the  signals  detected  by  seismic  analysts 
simply  cannot  be  retrieved  by  automatic  detector  algorithms 
based  solely  on  large  power  fluctuations.  Unfortunately, 
most,  if  not  all  conventional  automatic  seismic  signal  detec¬ 
tion  algorithms  fall  into  this  category.  For  this  reason, 
we  decided  that  simply  making  small  modifications  of  Unger's 
algorithm  would  not  result  in  a  sufficiently  improved  level 
of  signal  retrieval  performance.  To  achieve  that,  his 
decision  algorithm  would  need  to  be  generalized  to  operate 
on  the  basts  of  multivariate  amplitude  and  phase  measurements. 
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The  goal  of  utilizing  a  multivariate  automatic  signal 
editor  is  to  retrieve  signals  of  known  events  with  the  ef¬ 
fectiveness  of  a  seismic  analyst.  As  previously  discussed, 
magnitude  measurements  of  seismic  noise  are  highly  skewed. 
Noise  magnitudes  computed  as  the  logarithm  of  envelope  mea¬ 
surements  ,  as  is  done  by  Unger's  algorithm,  occur  most  fre¬ 
quently  at  the  RMS  level  of  the  seismic  noise.  The  standard 
deviation  about  the  most  probable  occurring  noise  magnitude 
is  0.5  magnitude  units  for  negative  deviations;  and  0.2  mag¬ 
nitude  units  for  positive  deviations. 

As  a  result  of  this  skewness,  small  seismic  signals  up 
to  0.5  magnitude  less  than  RMS  noise  can  be  expected  to  be 
larger  than  the  noise  in  10  to  20%  of  the  signal  windows 
examined.  This  provides  an  analyst,  skilled  in  recognizing 
signals  by  their  waveform  character,  a  reasonable  opportunity 
of  detecting  at  least  one  or  two  weak  signals  from  a  network. 
We  expect  that  by  discerning  change  of  signal  frequency  as 
well  as  amplitude  over  the  expected  duration  of  a  possible 
signal,  the  analyst  can  extract,  time,  and  measure  signals 
well  below  the  threshold  of  an  automatic  power  detector. 

Such  small  signals  would  always  be  missed  by  a  power  detec¬ 
tor  . 


Given  the  task  of  automatically  retrieving  signals  of 
known  events,  we  will  generalize  the  detector  so  it  can 
recognize  signals  based  on  a  multivariate  set  of  observa¬ 
tions  of  frequency  and  amplitude.  3y  doing  this,  there  is 
at  least  some  possibility  of  attaining  the  level  of  signal 
retrieval  performance  expected  from  a  seismic  analyst.  This 
is  expected  to  yield  as  much  as  0.5  magnitude  units  of  en¬ 
hanced  magnitude  capability  of  extracting  small  signals. 
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1.  Definition  of  a  Standard  Set  of  Multivariate  Seis¬ 
mic  Amplitude  and  Phase  Measurements 

The  basic  problem  of  defining  a  set  of  standard  multi¬ 
variate  measurements  is  to  encompass  the  physical  character¬ 
istics  of  earthquake  or  explosion  seismic  sources.  The  first 
step  in  achieving  this  is  to  transform  seismic  measurements 
to  a  multivariate  set  of  ground  motion  measurements.  In 
principal ,  this  theoretically  whitens  the  source  between  two 
source  dependent  frequencies.  By  sensing  the  range  of  peak 
frequencies  of  a  possible  signal,  an  optimum  filter  can  be 
automatically  designed  to  extract  the  signal. 

The  multivariate  set  of  ground  motion  measurement  will 
be  based  on  the  system  response  removal  function  now  being 
implemented  into  the  VS C  seismic  system.  The  seismic  record 
will  be  transformed  into  measurements  of  the  integral  of  the 
ground  displacement,  ground  displacement,  velocity,  accelera¬ 
tion,  and  the  derivative  of  the  acceleration.  From  Randall 
(1973),  the  asymptotic  behavior  of  earthquake  source  models 
in  terms  of  normalired  frequency,  x=f/f  ,  where  f  is  the 

w  w  * 

corner  frequency,  are  given  as  follows  for  the  displacement 
amplitude  spectrum  in  the  form  n(a>)  =Q(0)F  (x)  . 

a.  A  function  derived  by  Keilis-Borok  (1959) 
equated  to  the  low-  and  high-frequency 
asymptotes  intersecting  at  the  corner  frequency, 

F(x)  «  1  ,  (0<x<l) 

=  x'  2  ,  ( x>  1)  . 


r  T  r 
1  i  l 


n  - 

-  ± 


EXSCC,  INC. 


b.  The  function  considered  by  Brune  (1970), 


?(x) 


c.  The  spectral  shape  for  either  of  the  source 
models  of  Randall  (1966)  and  Archambeau  (1968), 


d.  A  function  given  by  Randall  (1973)  defined 
asymptotically  like  function  1  but  assuming  a 
w'J  high  frequency  behavior, 


F(x)  =  1  (0<xsl) 


e.  Two  models  given  by  Aki  (1967)  were  derived  by 
Haskell's  method  from  the  spatial  and  temporal 
correlation  of  the  velocity  of  a  fault  dis¬ 
location, 

F(x)  =  (1+ax2)  *  (l+x“)  * ,  w"  model 
F(x)  =  ( 1+ax")  (1+x2)  ,  wJ  model. 


r.  Mueller  (1969)  derived  the  displacement  ampli¬ 
tude  spectrum  of  explosion  sources  from  the 
Latter,  et  al .  (1959)  model, 

F ( x)  =  x1  (x6+a1x4+a,x2+a-) . 


andall  (19 "3)  showed  that  expressions  derived  to  compute 
eisnic  energy  and  characteristic  stress  are  independent  of 
ssumptions  as  to  source  model.  Thus,  all  of  the  models 


p  r  e  v  i o us iy  listed  o  r  o  v  i  d  e 
oroadband  ground  notion  fi 
characteristics  of  seismic 
mat ic ally  in.  Figure  I II -4. 


a  reasonable  basis  for  designing 
Iters  to  extract  the  spectral 
signals.  These  are  shown  sche- 


The  models  illustrated  in  Figure  MI-4  show  that  moments 
w_  tne  signal  spectrum  are  generally  whitened  over  a  portion 
-he  trequencv  band  controlled  by  the  source.  In  practice, 
_nia  et_ect  sometimes  will  be  masked  by  path  absorption  and 

.ecener  scattering  requiring  source  region-station  calibra¬ 
tion. 

t  h. 

ihe  k  spectral  moment  of  the  normalized  displacement 
amplitude  distribution  as  a  function  of  frequency,  F(x),  is 
given  by 

”  k 

4v  =  /  x  F (x)  dx  . 

0 


i ime  domain  estimates  of  can  be  obtained  by  taking  the  kL^ 
derivative  of  the  seismic  record,  detecting  and  timing  the 
signal,  and  measuring  the  magnitude  of  the  signal.  3y  mea- 
S— »j.ng  Jie  -requency  corresponding  to  envelope  peaks,  an 
estimate  can  be  obtained  of  the  frequency  band  of  the  k^ 
spectral  moment. 


nae 


requency 
should  in 
istics  of 


time 


and 

pri 

the 


domain  operators  which  whiten  the  signal  models 
-4  are  the  k  derivatives  where  k= - 1 , 0 , 1 , 2 , 3 . 
magnitude  measurements  of  these  derivatives 
nciple  completely  specify  the  source  character- 
seismic  event . 


:xc. 
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FIGURE  III -4 

ILLUSTRATION  OF  SEISMIC  SOURCE  MODELS 


The  set  of  k“n  derivatives  (-1,0 ,1,2, 5)  of  the  displace¬ 
ment  ground  motion,  each  reduced  to  measurements  of  magnitude 
anc  rrequencv ,  will  result  in  ten  multivariate  statistics 
representing  each  event.  To  utilire  these  measurements,  mag- 
n.-ude  and  rrequencv  measurements  of  each  derivative  need  to 
oe  calibrated  as  a  function  of  event  magnitude  since  the  cor¬ 
ner  frequencies  are  magnitude  dependent.  Also,  the  calibra¬ 
tion  is  needed  to  correct  for  path  absorption  and  site  charac¬ 
teristics  . 


dirricult  problem  is  anticipated  in  separating  the 
whitened  signal  bands  from  seismic  noise;  especially  for  the 
ground  displacement  (0L^  derivative)  and  the  integral  of  the 
g.ound  displacement  (-1  derivative) .  As  pointed  out  in 
-he  discussion  of  noise  models  in  Section  II,  the  existance 
or  nearly  periodic  peaks  mixed  with  nearly  white  Gaussian 
noise  will  tend  to  change  the  skewed  statistical  distribu¬ 
tion  or  noise  magnitudes  to  a  more  symmetrical  normal  dis¬ 
tribution  of  much  lower  standard  deviation.  If  the  observed 
rrequencv  of  noise  peaks  of  the  kr^  derivative  trace  are  out¬ 
side  the  anticipated  signal  frequency  band,  e.g.,  such  as 
m^cr oseisms  peaking  at  0.2  Hz.  Filtering  to  extract  the 
nearly  periodic  noise  component  and  applying  amplitude  and 
-hase  demodulation  techniques  to  minimize  the  interference 
or  microseisms,  it  should  be  possible  to  improve  measurements 
of  weak  signals. 


I  i  I  -  2  3 
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BASIC  GROUND  MOTION  MEASUREMENTS 

In  Section  II,  we  discussed  errors  associated  with  the 
use  oc  me  analytic  technique  o£  Unger  and  Fambach;  espe¬ 
cially  the  problem  of  making  precise  estimates  of  the  fre¬ 
quent-  01  signals .  In  that  case,  the  error  increases  without 
-..mi-  as  -he  envelope  approaches  zero.  More  generally,  if 
envelope  changes  over  the  time  interval  of  the  frequency 
es  t.ma.e ,  -hen  the  rrequencv  estimate  is  biased  or  will  lack 
rrecision . 


.-or  this  reason,  Unger's  po int -by -point  technique  of 
measuring  signals  as  the  logarithm  of  the  envelope  contained 
in  a  4  second  window  is  modified  as  follows . 

•  Seismic  records  are  transformed  into  sequences  of 
measureable  envelope  peaks 

•  Data  associated  with  each  envelope  peak  are 

magnitude  of  ground  motion 

frequency 

arrival  time. 

Ir  a  seismic  peak  is  encountered,  we  will  have  timed 
the  envelope  minimum  preceding  the  peak  and  that  following 
the  peak.  3y  symmetry  considerations,  it  can  be  shown  from 
Rice  (1954)  that  most  probable  occurrences  of  envelope  mini¬ 
ma  or  a  Gaussian  time  series  are  approximately  0.5  of  the 

RMS  of  the  time  series;  of  envelope  maxima,  1.5  times  the 
RMS  . 


To  distinguish  signal  envelope  peaks,  the  envelope  level 
o:  tne  maxima  must  be  1.5  to  5  times  the  adjacent  minimum 
values .  Arter  applying  this  criteria,  the  magnitude  of  the 
seismic  signal  peaks  will  be  computed  as  the  common  logarithm 
p-us  a  transmission  3-factor  to  correct  for  propagation  from 
a  source  of  known  location. 

.ne  problem  or  determining  accurate  frequency  measure¬ 
ments  is  minimized  by  avoiding  determinations  at  times  where 
there  are  large  negative  magnitude  fluctuations.  The  fre¬ 
quency  measurements  are  optimized  by  determining  the  fre¬ 
quency  at  maxima  of  the  envelope  peaks.  Figure  III-5  shows 
an  envelope  peak  occurring  between  the  ith  and  (i+l)th  data 
point  on  a  seismic  record. 

3v  interpolation,  the  envelope  peak  shown  in  Figure 
III-5  can  be  modeled  by  a  cosine  modulation. 

E(t)  =  EQcoswit  . 

Tms  adequately  represents  the  time  variation  of  the  envelope 
m  tne  immediate  neighborhood  of  the  interpolated  maximum  at 
point,  x,  shown  in  Figure  III -3.  Similarly  by  interpolation, 
tne  seismic  data,  x(t) ,  is  interpolated  between  points  i  and 

i-  i ; 

X(t)  s  E( t)  sin  (wQt  *  $  ) 

=  :  cos  w  t  [cos  b  s  in  w  t  *  sin  6  cos  w  t] 

OX.  0  0  0  0  J 
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puK  i  un 


=  *  =  0  {cos  [sin(wQ+wz)t  *  s  in(wQ -w£)  t] 

-  sin  o Q  [cosCwQ-wA)t  +  cos (w-w£) t  ]  } . 

.\e  see  that  even  under  these  ideal  conditions  of  measuring 
rrequency  at  the  envelope  peak,  time  varying  phase  angle 
estimates  of  the  frequency  could  be  seriously  in  error  due 
to  envelope  modulations  of  angular  frequency,  w^.  Nonethe- 
_ess ,  at  the  interpolated  maximum  envelope  peak  an  unbiased 
rrequency  estimate  can  be  obtained  by  combining  envelope  and 
phase  angle  measurements  with  the  time  derivative  of  the  seis¬ 
mic  trace,  as 


x'  Co) 


dxCt)  | 
at  1 t=0 


w  E  cos <j>  . 

oo  o 


Since  the  Hilbert  transform  or  quadrature  filtered  seismic 
trace,  y(t),  is  x(t)  phase  shifted  90°;  the  time  derivative 
is 


y 1  Co)  = 


dy(t)  | 

at  1 1=0 


•w  E  s in® 
oo  o 


..ius,  an  unbiased  estimate  of  the  frequency  is  obtained, 


C 1 1 1  -  4 ) 


at  the  interpolated  maximum  envelope  peak.  This  requires 
--me  derivatives  ot  the  seismic  trace  and  its  Hilbert 
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rods 


■  srcrm  and  measurements  of  the  envelope  at  that  point, 
.pplymg  the  preceding  analysis,  it  can  be  shown  that  by 
ling  envelope  modulations  at  inflection  points  as  a  con- 
t  envelope  value  EQ  plus  a  sinusoidal  modulation 
nw,t ,  the  slope  at  such  a  point  is  given  as 


x'(o) 


E  w  coso 

O  0 


AE  w . s ino 
ox.  o 


omce  AEqw  is  the  derivative  of  the  envelope  function  at  t=0 , 
we  see  that  computing  frequencies  at  points  other  than  the 
peak  or  the  envelope  function  yields  bias  due  to  the  slope 
of  the  envelope  function.  Thus,  equation  (III-4)  should  be 
used  to  compute  frequency  at  interpolated  envelope  peaks  to 
avoid  bias  due  to  envelope  modulation. 

The  problem  of  phase  modulation  caused  by  random  inter¬ 
ference  of  other  frequencies  is  another  source  of  error. 

Near  peaks  and  troughs  of  the  random  phase  angle  modulation, 
the  correct  frequency  will  be  given  by  equation  (III-4). 

Near  rero  crossings,  maximum  positive  or  negative  deviations 
rrom  the  correct  frequency  are  caused  by  phase  modulations. 
Unrortunately ,  if  we  only  measure  frequency  at  envelope  peaks  , 
we  are  exposed  to  this  source  of  error.  3ut ,  if  we  initially 
estimate  the  frequency  by  equation  (III-4),  data  conditioning 
can  be  carried  out  to  minimize  this  source  of  error. 

The  time  series  containing  the  envelope  peak,  can  be 
interpolated  and  filtered  by  a  bandpass  filter  of  sufficient 
bandwidth  to  pass  the  peak  without  distortion.  This  is  given 
by  Rice  (1954)  both  for  an  ideal  (square)  bandpass  filter  and 
an  equivalent  bandwidth  Gaussian  bandpass  filter.  Such 
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: i- -ermg  about  tne  initially  estimated  frequency  will  re- 
auce  phase  modulation  errors. 

rurther  reduction  of  phase  modulation  errors  is  possible 
oy  averaging  the  frequency  measurements  over  a  time  span 
covering  that  portion  or  the  time  series  exceeding  the  most 
probable  occurring  envelope  from  some  hypothetical  Gaussian 
input  giving  rise  to  the  observed  envelope  peak.  This  can 
-e  aetermined  approximately  and  robustly  by  measuring  the 
frequency  at  points  about  the  envelope  peak  which  are  above 
“he  median  envelope  level .  To  avoid  errors  due  to  envelope 
modulation,  the  filtered  data  x(t)  and  the  Hilbert  transform 
■"  C c )  are  uransrormed  to  minimize  the  effect  of  envelooe 
modulation . 

uCt)  =  xCt) CEQ/E  Ct) ) 

vCt)  3  y Ct) (E0/E(t))  . 

The  time  series  u(t)  and  v(t)  are  substituted  for  x(t)  and 
y(t)  in  equation  (111-4) .  Since  envelope  modulation  effects 
are  reduced  by  this  process,  the  frequency  can  be  measured 
at  all  points  m  the  neighborhood  of  the  envelope  peak  above 
“he  median  of  the  envelope  function  defining  the  peak.  By 
applying  this  procedure  to  measure  the  frequency  associated 
with  each  observed  peak  envelope  we  expect  to  obtain  more 
precise  estimates  of  frequency  than  those  previously  obtained 
°v  Unger  and  others.  Our  frequency  estimates  will  be  more 
precise  and  minimally  biased  by  minimizing  errors  caused  by 
envelope  and  phase  modulation. 


HI 
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By  measuring  seismic  ground  motion  with  the  above  tech¬ 
nique,  we  expect  to  obtain  accurate  magnitudes  corresponding 
to  observed  envelope  peaks.  In  addition,  estimates  of  the 
dominant  frequency  of  ground  motion  will  be  associated  with 
each  observed  envelope  peak.  The  arrival  time  of  each  peak 
will  be  taken  as  the  time  of  the  first  envelope  minimum 
preceding  the  observed  envelope  peak. 

The  input  to  such  a  basic  ground  motion  measure  package 
would  be  the  n  ‘  derivative  of  ground  displacement 
(n= - 1 , 0 , 1 , 2 , 3) ;  the  output  of  a  sequence  of  envelope  peaks, 
each  specified  by 

•  Magnitude 

•  Frequency 

•  Arrival  time. 

Seismic  records  will  be  converted  sequentially  into  five 
derivative  ground  motion  time  series.  In  turn,  each  of  these 
derivative  ground  motion  traces  will  be  transformed  into  a 
corresponding  set  of  the  basic  ground  motion  measurements 


J  . 


.N'OISE^DATA^CONDITIONINGGROt^^D  M°TI°N'  MEASUREMENTS  °E 


The  first  step  in  associating  basic  ground  motion  mea¬ 
surements  with  signals  from  a  known  event  is  to  distinguish 
possible  signals  rrom  magnitude  peaks  associated  with  seis¬ 
mic  noise,  each  record  accessed  after  locating  and  timing 
a^e.smic  event  will  contain  approximately  one  or  two  minutes 
or  seism: c  noise  data  starting  about  2  or  3  minutes  before 
tne  expected  signal  first  arrival  time  and  ending  about  1 
minute  oerore  expected  signal  arrival  time.  The  purpose  of 
applying  a  1  minute  guard  gate  is  to  avoid  mixing  signal  data 
*-.h  that  wmch  is  to  be  interpreted  as  noise.  Since  ordered 
statistics  will  be  used  to  derive  the  noise  parameters  needed 
-o  detect  signals,  mixing  of  a  small  fraction  of  signal  in¬ 
tonation  with  noise  should  not  seriously  affect  the  deter¬ 
mination  of  noise  parameters.  The  advantage  of  the  ordered 
statistics  approach  is  their  robustness.  3y  neglecting 
extreme  large  and  small  deviations  effects  produced  bv  inter¬ 
mittent  malfunctions  such  as  spikes,  signals  mixed  in' noise, 
s-c.,  will  be  minimized. 

.-allowing  the  Gaussian  noise  model  described  in  Subsec¬ 
tion  A,  tne  magnitude  level  bounding  205  of  the  smallest  mag¬ 
nitudes,  o65  and  38%  are  used  to  derive  the  most  probable 
occurring  noise  magnitude,  as  well  as  the  standard  deviation 
=i  positive  and  negative  magnitude  deviations.  Under  the 
ideal  Guassian  assumption,  the  negative  standard  deviation  is 
expected  to  be  about  2.5  times  larger  than  the  positive  devia¬ 
tion  and  the  most  probable  envelope  peak  magnitude  is  expected 
-o  occur  at  1.3  times  the  RMS  of  the  noise.  If,  however', 
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•  e.a^i.eiy  oroadband  Gaussian  noise  is  mixed  with  an  extremely 
narrows  ana  (almost  periodic)  component  (Rice,  19  54),  the  ratio 
c.  negative  to  positive  standard  deviation  is  close  to  one. 
ohoulc  this  condition  De  detected,  and  should  a  narrowband  of 
rrequencies  oe  indicated  by  the  ordering  of  frequency  star.is- 
-ius  (e.g.,  ,0?j  or  the  observed  frequencies  between  the  octave 

"nu  ^  and  0.4  H«) ,  then  the  data  must  be  conditioned  to  re¬ 
move  this  noise  source. 


^ata  conditioning  applied  to  demodulate  useful  seismic 
cata  rrom  data  containing  an  interfering  nearly  periodic 
noise  component  is  illustrated  by  Figure  III-6.  The  equations 
ror  separating  a  desired  broadband  process,  Ax(t),  from  the 
mixed  data  xy(t)  are  as  follows. 


3y  insPection  of  the  illustration  and  the  definition  of 
variables  in  Figure  III-6,  we  have  the  following  equations 
for  demodulation  of  the  small  broadband  data  component  Ax(t) . 


o'(t)  =  TT-y(t)  -  er(t) 


By  law  of  cosines, 

AR(t)  =  [Ru(t)  +  Rf(t)  -  :.Ry(t)  Rf( t)  Cos  (A0)J ^ 

where 


ie(t)  =  s  ft)  -  e-(t) 

H  £ 


AR(t) “  +  Rr(t)"  -  R  (t)  “ 

cosv(t)  =  - — _ _  ^ 

-^■-1")  AR(t) 


\ 

I 


AyCt) 


Similarly  (x^,  y.,  R^,  and  0^)  are  complex  seismic  data 
points  rrom  the  filtered  narrowband  process 

envelope  of  demodulated  data  (effect  of  interfering 
narrowband  component  removed) 

3  Phase  angle  of  demodulated  data 

-x(c)  Demodulated  broadband  data  obtained  bv  removing  the 
strong  interfering  spectral  component  xc(t)  from 

untiltered  data  x  (t) 

u^  J 


FIGURE  1 1 1  -  6 

SEPARATION  OF  A  BROADBAND  PROCESS  FROM  THE  EFFECT  OF  A 
STRONG  INTERFERING  NEARLY  PERIODIC  COMPONENT 
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•r*e  cesired  demodulated  broadband  data  are  aiven  as 


-xCt)  =  AR(t)  coscCt) 


CHI  -  5) 


t_mp_>  get-ing  rid  ot  the  strong  almost  periodic  component  by 
**“Sa*  -altering  is  orten  insumcient  because  the  modulation 
o.  v.he  small  broadband  noise  component  effectively  broadens 
tr.e  tandvidth  of  the  interfering  narrowband  process.  This 
ma.\es  it  difricult  to  see  signals  at  frequencies  near  the 
...c5..er*ng  spectral  peak.  This  can  be  seen  intuitivelv  by 
examining  rigure  I II -6.  The  broadband  data  containing  much 
sher  rrequencies  moves  around  point  P  many  cycles  while  the 
strong  low  frequency  interference  component  moves  point  P 
very  slowly.  This  has  the  effect  of  modulating  the  envelope 
Rc ,  between  R£-dR  and  R£+AR.  It  also  modulates  the  phase  ,0, 
oetween  G^itan*  1(AR/'R£)  .  These  broaden  the  bandwidth  of  the 
spectral  peak.  A  broader  band  filter  is  required  to  remove 
this  periodic  component. 


Having  thus  conditioned  the  data,  we  are  in  a  position 
of  generating  and  interpreting  noise  statistics.  At  least  as 
a  reasonable  approximation,  the  noise  can  be  interpreted  as  a 
Gaussian  process. 


On 

termite 
is  done 


that  basis,  we  use  ordered  statistics  to  robustly  de- 
uhe  most  probable  occurring  noise  magnitudes.  This 
separately  tor  each  of  the  n^  derivative  seismic 
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A 1th ough 
normalic  at 


the  distributions  are  skewed,  we  use 
on  which  divides  positive  deviations 


a  method 
by  a 
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positive  standard  deviation; 
star.aard  deviation.  For  the 
this  ultimately  reduces  all 
to  a  homogeneous  set  of  unit 


negative  deviations  by  a  negative 
purpose  of  multivariate  analysis, 
of  the  signal  detection  variables 
normal  noise  statistics. 


.-or  each  n  derivative  ground  displacement  seismogram, 

’,e  5  ^a* »  wi  uh  a  set  or  three  Basic  Ground  Motion  measurements, 
.r.ese  are  the  arrival  time  of  an  envelope  peak,  the  magnitude, 
o_  -he  envelope  maxima,  and  the  dominant  frequency. 


rrom  ordered  statistical  analysis  of  the  noise  data 
preceding  the  window  containing  signals,  we  generate  statis¬ 
tics  or  magnitude  and  frequency  fluctuations . 


One  of  the  main  purposes  of  using  a  multivariate  detector 
to  time  and  measure  signals  is  to  improve  our  ability  to  auto¬ 
matically  time  and  measure  small  signals.  These  presumably 
could  be  detected  as  statistically  significant  deviations  of 
uhe  frequency  or  magnitude  from  values  expected  of  seismic 
noise.  These  would  be  expected  to  persist  for  a  duration  of 
0.5  to  10  seconds  (the  expected  duration  of  a  seismic  phase)  . 

For  that  purpose,  we  separately  determine  relationships 
_or  positive  and  negative  deviations  from  the  most  probable 
occurring  noise  magnitude.  Let  AM  and  f ^  be  the  normalised 
magn-wude  fluctuation  and  corresponding  frequency  of  the  i“^ 
observed  noise  envelope  peak.  Equation  (III-6)  shows  the 
equations  for  determining  a  statistical  noise  prediction 
model  for  positive  magnitude  fluctuations. 
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(III-6) 


where  1  is  a  vector  of  one's; 
tions;  and  £_+,  the  associated 
deviations . 


positive  magnitude  devia- 
frequency  of  positive  magnitude 


The  purpose  of  (III-6)  is  to  derive  a  systematic  rela¬ 
tionship  between  positive  noise  magnitude  deviations  and 
their  associated  frequency.  The  relationship  is  used  to  pre 
diet  the  most  probable  magnitude  of  noise,  m,  a  constant  vec 
tor  .  Its  value  is  equal  to  the  most  probable  noise  magnitude 
aerived  from  ordered  statistical  analysis  of  the  observed 
noise  magnitudes. 

Similarly  equation  CHI-7) ,  determines  a  noise  predic¬ 
tion  error  operator  b  for  negative  magnitude  fluctuations 
vmcn  utilizes  the  frequency  of  observed  noise  peaks  to 

minimize  their  deviation  from  the  most  probable  observed 
no ise . 


1  1  1 

b. 

1 

1_  AM  ~  1  o  g  f 

0  2 

b- 

_  j 

= 

m 

i  1  T 

+ 
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The  physical  rationale  for  this  approach  is  seen  in  Figure 
III-".  Unger  (19"")  observed  a  linear  relationship  between 
dominant  frequency  and  phase;  or  equivalently,  the  frequency 
standard  deviation  or  bandwidth  of  seismic  signals.  This  is 
shown  schematically  in  the  hypothetical  relationship  in 
Figure  1 1 1 -  7 .  Note  that  negative  magnitude  fluctuations  are 
postulated  as  broader  bandwidth  (higher  variance)  pulse 
spectra.  The  possible  association  of  higher  magnitude  fluc¬ 
tuations  as  more  coherent  lower  frequency  noise  pulses  of 
equivalent  seismic  energy  seems  to  be  a  plausible  concept 
for  the  noise  prediction  error  model  of  equations  ( 1 1 1  -  6 ) 
and  (III  -  7)  . 

Given  that  the  noise  prediction  error  model  works  as 
expected,  it  will  significantly  reduce  the  variance  of  noise 
magnitude  fluctuations.  This  results  in  more  reliable  extrac¬ 
tion  of  signals  of  known  seismic  events.  Even  small  magnitude 
fluctuations  caused  by  weak  signals  can  then  be  more  reliably 
associated  with  known  events,  timed,  and  measured.  The  ap¬ 
plication  of  equations  ( 1 1 1  -  6 )  and  (III-7)  is  expected  to  re¬ 
sult  in  sufficient  reduction  of  the  standard  deviation  in  the 
magnitude  of  noise  peaks,  that  small  signal  pulses  greater 
than  two  standard  deviations  of  the  noise  peaks  will  reliably 
be  detected  by  Unger's  algorithm  discussed  in  Section  II. 

An  important  part  of  our  signal  extraction  strategy  is 
tc  use  frequency  measurements  as  a  dependent  variable  for 
detecting  weak  signals  which  are  as  much  as  a  half  magnitude 
below  the  noise  level.  We  do  this  by  applying  ordered  sta¬ 
tistics  to  measurements  of  the  frequency  of  noise  peaks. 
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3  Most  frequently  observed  magnitude  of  noise 
3  Noise  magnitude 
3  M-m  if  Sgn  (M-m)  >  0 
3  M-m  if  Sgn  (M-m)  <  0 

3  Log  of  dominant  frequency  of  envelope  peaks 


FIGURE  III- 7 

YPOTHETICAL  MODEL  SHOWING  MAGNITUDE  FLUCTUATIONS  OF  NOISE: 
HIGHER  FREQUENCY  ENVELOPE  PEAKS  ASSOCIATED  WITH  BROADER 
BANDWIDTH  AND  LOWER  AMPLITUDE;  LOWER,  WITH  NARROW 
3ANDWIDTH  AND  HIGHER  AMPLITUDE 


m 

M 
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AM- 
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Per  that  purpose,  we  separately  determine  relationships 
for  positive  and  negative  deviations  of  the  frequency  of 
noise  peaks .  Since  we  do  not  have  a  statistical  model  of 
frequency  fluctuations  to  draw  on,  we  assume  that  the  dis¬ 
tribution  is  skewed.  The  median  frequency,  F,  is  taken  as 
the  most  probable  occurring  frequency  of  noise  envelope  peaks 
Positive  and  negative  standard  deviations  are  then  at  the 
and  1 5 wl  percentile  of  observed  frequencies  of  noise 
peaks  ordered  from  lowest  to  highest  frequency.  Let  If.  be 
the  observed  frequency  fluctuation  and  NL  the  associated 
magnitude.  Following  the  methodology  resulting  in  equation 
(III-6),  we  have  the  corresponding  equations  for  observed 
positive  frequency  fluctuations. 


1  If*  NT 


( 1 1 1  -  8) 


Allowing  the  methodology  of  equation  (III-7),  we  have  a  cor- 
•esponding  equation  for  observed  negative  frequency  fluctua- 


:rcns  . 


1  a  - 


C 1 1 1  -  9 ) 


1 1 1  -  4  1 


I 
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ihe  purpose  or  (III-S)  and  (III-9)  is  to  derive  a  sys¬ 
tematic  relationship  between  frequency  fluctuations  and  the 
magnitude  measurement  associated  with  each  observed  frequency 
r-uctuation.  The  physical  motivation  was  shown  in  Figure 
HI-'.  This  procedure  is  expected  to  minimice  the  variance 
of  frequency  fluctuations  of  noise  peaks.  This  technique 
makes  it  possible  to  detect  weaker  signals. 

For  some  signals  we  expect  that  it  will  be  preferable 
to  time,  extract,  and  measure  signals  based  on  significant 
magnitude  rluctuations .  In  other  cases,  we  expect  it  will 
be  perrerable  to  use  frequency  measurements  as  the  dependent 
variable  ror  this  purpose.  Unger  (1973)  showed  that  phase 
mgle  criteria  is  at  least  6  dB  more  sensitive  in  detecting 
’■veax  signals  than  is  amplitude.  Since  frequency  is  a  Dhase 
modulation  measurement,  it  should  be  a  more  sensitive  cri¬ 
teria  ror  detecting  those  signals  of  significantly  different 
rrequency  than  is  noise.  For  example,  it  is  commonly  ob¬ 
served  that  regional  Lg  phases  can  be  most  easily  recogniced 
by  abrupt  rrequency  changes.  For  these  reasons,  we  will  ap¬ 
ply  both  tne  magnitude  and  the  frequency  criteria  to  time, 
extract,  and  measure  weak  seismic  signals. 

Tne  operators  a,  b,  c,  and  d,  in  equations  (III-6) 
through  (III-9)  are  derived  from  the  analysis  of  noise  data 
provided  at  the  front  end  of  each  record.  These  are  optimum 
operators  ror  removing  the  effects  of  frequency  fluctuations 
from  magnitude  measurements  and  magnitude  fluctuations  from 
frequency  measurements.  Positive  and  negative  fluctuations 
are  separately  treated.  Letting  the  matrix  operator  G  or  H, 
represent  noise  observations  and  the  vectors  £  or  h,  one  of 
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the  desired  operators  for  removing  the  effect  of  fluctuations; 
we  can  determine  a  least  squares  solution  as  follows.  For 
estimating  magnitudes, 

l  =  (GTG) _1  GT  m  ,  (III-10) 

for  estimating  frequency, 

h  =  (HTH) _1  HT  F  .  C I II - 11) 

As  was  done  for  scaling  discriminants  in  the  Event 
Identification  Experiment,  more  robust  determination  of  £  and 
h  can  be  obtained  from  small  samples  subject  to  occasional 
large  errors.  This  is  done  simply  by  altering  the  rules  for 
vector  dot  multiplication  by  replacing  x*£  with  N  times  the 
median  of  {x^y^},  where  i*l,2,...,N.  Since  matrix  multiplica¬ 
tion  involves  a  set  of  vector  dot  products,  robust  determina¬ 
tions  involve  a  set  of  sort  operations  for  median  determina¬ 
tions.  It  is  noted  that  if  statistical  deviations  are  normal¬ 
ly  distributed  and  N  is  very  large,  robust  determinations  of 
the  operators  are  equivalent  to  least  squares  estimates. 

The  purpose  of  applying  the  operators  to  measurements  of 
envelope  peaks  is  to  obtain  minimum  noise  variance  estimates 
of  magnitude,  m,  and  frequency,  F.  These,  obtained  from  each 
of  the  n  4  derivative  ground  motion  measurements ,  are  desig¬ 
nated  m  ,  F  .  The  weak  signal  extraction  problem  gets  down 
to  distinguishing  the  ten  element  multivariate  vector 
(m_,,  F_1,  m  ,  FQ,...,m_,  F-) .  It  is  perhaps  worth  noting 
that  m-  versus  F-  has  been  used  as  an  effective  discriminant 
(referred  to  as  the  third  moment  method)  between  earthquakes 


and  explosions.  I 
multivariate  signa 
types  of  events  as 


t  is  not  unreasonable  to  expect  this 
1  extraction  process  to  distinguish  various 
well  as  distinguish  signals  from  noise. 


E.  MEASUREMENTS  ON  ENVELOPE  PEAKS;  POST  -  FILTERING 


i.  Statistically  Independent  Estimates 

In  Subsection  D,  we  derived  correction  operators  g  and 
n  ror  obtaining  m  and  F.  m  is  the  magnitude  fluctuation  as¬ 
sociated  with  an  envelope  peak  corrected  for  its  apparent 
zrequency  dependence;  F  the  frequency  fluctuation  corrected 
ror  its  magnitude  dependence.  The  corrections  were  designed 
'■o  minimize  the  variability  of  noise  fluctuations.  The  pur¬ 
pose  of  this  correction  was  to  maximize  the  sensitivity  for 
.automatically  detecting,  timing,  and  measuring  weak  signals 
occurring  in  some  unknown  position  of  the  accessed  time 
windows  of  a  known  event.  Any  statistically  significant  de- 
^  *at.ons  of  m  and  F  from  values  expected  for  noise  can  then 
be  interpreted  as  possible  onsets  of  seismic  signals. 

The  rirst  step  in  analyzing  noise  measurements  of  m  and 
r  1S  t0  normalize  m  and  F  to  unit  normal  variates  by  applica¬ 
tion  or  ordered  statistics.  3y  subtracting  the  most  probable 
occurrence  rrom  each  variate  and  dividing  by  the  appropriate 
positive  or  negative  standard  deviation,  m  and  F  are  reduced 
“ °  a?Pr°ximate  unit  normal  statistical  variates  Z  (m)  and  Z[F)  . 
ine  observations  of  Z(m)  and  Z(F)  are  operated  upon  by  rota¬ 
tional  trans rormation  operator,  Q,  which  transforms  them  to 


.  and  Z'(Fj  which  are  statistically  independent.  This 

operation  is  shown  by  equation  (III-12). 


where  Q  is  a  2x2  matrix.  By  taking  the  transpose  of  the  ob¬ 
servation  matrix  and  applying  it  to  both  sides  of  equation 
(III-12)  we  obtain  a  least  squares  determination  of  Q  . 


var{Z(m)}  cov{ Z(m)  , Z (F) } 

RQ  =  Q 

cov{ Z (m) ,  Z(F) }  var(ZCF)} 


1.0  0 

Cm- is) 

o  1.0 


The  solution  for  Q  is  simply  obtained  as  the  inverse  of 
the  covariance  matrix,  R,  of  magnitude  and  frequency  fluctua¬ 
tions.  Q,  by  definition,  and  when  applied  to  observations 
Z(m)  and  Z(F) ,  produces  Z ' (m)  and  Z'(F)  which  are  statistically 
independent  variables  for  detecting  signals  by  significant 
magnitude  or  frequency  fluctuations  from  noise.  Since  m  and 
F  have  been  pre-conditioned  to  remove  any  obvious  dependencies, 
it  is  expected  that  R  ^  will  be  a  very  stable  operator.  The 
rotation  operator  Q  is  applied  to  remove  any  apparent  depen¬ 
dency  between  m  and  F  so  that  the  assumption  of  independent 
Gaussian  statistics  is  val:.d. 


As  desired,  we  now  have  formulated  two  independent 
.'uawistics  tor  detection  signals.  I'(m)  is  an  approximate 
unit  normal  statistic  for  detecting  signals  by  their  unusually 
large  magnitude  fluctuations  from  noise.  Z'(F)  detects  small 
s-gnals  oased  on  unusually  large  changes  in  frequency  from  the 
most  probable  rrequency  occurrence  of  noise. 

2.  Bandwidth  of  Envelope  Peaks  and  Post-Filtering  of 
Signals 

Cur  rationale  in  separately  searching  for  signals  in  n^ 
derivative  ground  motion  traces  of  a  record  containing  signals 
rrcm  a  known  event  was  illustrated  schematically  by  Figure 
iiI-4.  The  basic  idea  was  to  whiten  signals  between  pairs  of 
corner  frequencies  governing  the  spectrum  of  a  signal.  For 
example,  a  3rune  model  source  can  be  considered  to  be  a  low 
bandpass  white  source  of  ground  displacement  (0^  derivative) 
up  to  near  the  comer  frequency  and  a  high  bandpass  white 
source  of  ground  acceleration  (2nc^  derivative)  above  the 
corner  frequency.  Suppose  such  a  source  was  measured  with 
seismic  data  which  are  band  limited  between  f^  and  f7.  Then 
the  bandwidth  of  ground  displacement  would  be  expected  to  be 
**“  (cepth)  =  fc"f^»  where  f^  is  the  corner  frequency;  of 
s*ound  acceleration,  Af(accel)  =  By  separately  analyc - 

mg  all  five  motion  derivatives  between  -1  and  3  we  expect  to 
ob  mm  physically  meaningful  band-limited  white  representation 
cr  the  source  in  different  portions  of  the  spectral  band.  3v 
measuring  one  or  more  detected  envelope  peaks;  their  magnitude, 
rrequency  and  bandwidth;  we  expect  to  obtain  useful  source- 
related  information. 
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at  which  e ravel 
Guassian  input 
cf  noise  is  At 
noise  process , 


•0  analyced  the  problem  of  determining  the  rate 
ope  peaks  are  passed  by  a  band  limited  white 
If  the  average  time  between  envelope  peaks 
seconds,  then  the  effective  bandwidth  of  the 
Af,  is  derived  by  Rice  as 


—  U.  -  1.30  .  (.111-14) 

.nus,  ic  we  observe  noise  envelope  peaks  occurring  on  an  aver¬ 
age  interval  of  Atn  seconds,  then  the  effective  bandwidth  of 
tne  noise  process  can  be  approximated  as  1.56/At  .  Since 
^nger  (19/8)  demonstrated  that  seismic  noise  envelope  measure¬ 
ments  can  be  reasonably  modeled  as  a  Gaussian  process,  this 

should  be  a  reasonable  assumption  for  samples  of  noise  of  long 
duration. 

A  more  interesting  parameter,  which  can  be  applied  to 
observations  of  individual  envelope  peaks  is  the  relationship 
between  the  measured  pulse  width  of  envelope  peaks.  Our  basic 
envelope  peak  measurements  include  the  start  time  as  the  first 
significant  minimum  preceding  the  envelope  peak  and  the  end 
“ime  aa  the  first  significant  minimum  following  the  envelope 
peak.  Defining  the  measured  pulse  width  of  noise,  T  ,  as  the 
-an -.age  pulse  width  of  noise  envelope  peaks,  then  we  can  de¬ 
termine  a  constant,  a,  which  relates  measured  pulse  width  to 
A,nas  a.  .  This  yields  a  relationship  between  effective 
bandwidth  and  measured  pulse  width  as 


(1 1 1 - 15) 


The  purpose  of  measuring  the  effective  bandwidth  of  an  envelope 
peak  or  sequence  of  envelope  peaks  is  to  post-filter  ^eak  sig¬ 
nals.  The  filter  will  enhance  the  signal-to-noise  ratio  of 
weak  signals  as  part  of  the  signal  extraction  process  and  to 
extend  and  to  enhance  the  quality  of  subsequent  signal  mea¬ 
surements.  Figure  1 1 1 - S  shows  examples  of  post-filters  ap¬ 
plied  to  extracting  weak  signals.  Af  is  the  full  bandwidth  of 
a  filter  designed  to  pass  weak  signals.  For  the  trapezoidal 
filter,  Af,  is  obtained  from  equation  (III-14).  For  the  Gaus¬ 
sian  filter,  f,  it  was  derived  by  Rice  C.  19  54)  and  is  shown 
by  equation  CHI-16)  . 

AfA-r  =  0.67  .  CIII-16) 

Rice's  derivation  was  on  the  basis  of  the  bandwidth  passing 
the  same  average  amount  of  power  as  an  ideal  bandpass  filter. 

The  approach  toward  applying  post- filtering  is  to  initially 
set  a  high  false  alarm  rate  threshold  to  detect  weak  signals. 

If  such  weak  signals  are  detected,  post-filtering  will  be  ap¬ 
plied.  This  will  be  followed  by  applying  a  more  exacting 
threshold  permitting  acceptable  network  operating  characteris¬ 
tics  for  extracting  signals. 


IGURH  III 


DETECTING  COMPLEX  SIGNALS 


1.  Random  Signal  Characteristics 

Due  to  strength  heterogeneity  of  the  source  region  and 
other  complicating  factors,  the  signal  does  not  always  prop¬ 
agate  as  a  single  envelope  peak.  Such  a  source  can  be  more 
generally  characterized  as  a  randomly  distributed  source 
which  can  be  modeled  as  a  cluster  of  delayed  transmissions 
over  the  signal  duration,  D.  Over  a  longer  time  scale  the 
source  density  decays  and  the  signal  gradually  fades  out. 

Observations  by  Aki  (1969)  of  long-time  terms  of  coda 
well  below  the  signal  level  indicate  that  such  coda  can  be 
explained  by  source  and  receiver  back-scattering  at  frequen¬ 
cies  dependent  only  on  travel-time  from  the  source.  This 
effect  can  be  viewed  as  an  illumination  of  the  medium  near 
the  earth's  surface,  all  along  the  path  from  the  seismic 
source,  and  due  to  multiple  scattering  of  waves  observed  well 
after  the  signal  arrival  and  at  much  lower  levels  than  the 
signal.  The  time-frequency  scattering  function  observed  near 
the  source  by  Aki  indicates  that  the  effect  fades  rapidly  due 
to  spreading,  dispersion,  and  absorpt  on.  There  he  observed 
a  fading  amplitude  with  the  dominant  period  shifting  rapidly 
to  longer  periods  with  time  elapsed  from  the  event  origin 
time.  The  same  dominant  periods  are  seen  at  a  simultaneous 
time  delay  from  the  origin  time  of  the  source,  at  different 
receiver  sites,  at  varying  distance  from  the  source. 
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.n  completely  cirrerent  coda  model  is  relevant  to  the 
cetection  or  signals.  It  stems  from  the  concept  of  a  seismic 
source  distributed  randomly  in  time  and  space.  Here  our  con¬ 
cern  is  with  determining  a  B-factor  correction  applicable  to 
fading  coda  measurements.  IV e  view  the  signal  as  a  uniformly 
cistribured  random  source  modeled  as  shot  noise  over  the  sig¬ 
nal  duration,  D.  This  will  be  characterired  as  a  cluster  or 
envelope  peaks  of  approximately  the  same  magnitude  and  of 
similar  spectral  characteristics.  For  times  greater  than  D, 
we  continue  to  see  delayed  replicas  of  the  signal,  the  only 
cirrerence  being  that  the  magnitude  decreases  after  the  sig¬ 
nal  duration,  D,  has  elapsed  and  then  fades  gradually  with 
time.  This  differs  from  Aki ' s  model  in  that  it  pertains  to 
snort-term  time  delays  at  initial  magnitudes  close  to  that 
or  tne  signal.  This  particular  model  was  used  to  simulate 
earthquakes  by  Shoup  and  Sax  (1974).  It  realistically  pre- 
cicted  the  time-span  of  visible  signals  over  a  wide  range  of 
magnitudes.  The  fading  shot  noise  coda  model  B-factor  cor¬ 
rection  is  given  by  equation  (III-17). 


d3c(t)  =  a  log  (t/D)  +  c, 


(t>d) 


(I  II  -  1 7) 


where  t  is  the  elapsed  time  from  the  start  of  the  signal, 
.ms  is  in  contrast  tr  Aki's  model,  where  he  computes  decay 
with  time  as  the  elapsed  time  from  the  origin  time  of  the 
event.  This  model  is  considered  valid  only  for  short-term 
cecay  or  signals  where  Aki's  model  is  valid  for  the  long¬ 
term  decay  occurring  well  after  the  arrival  of  the  signal. 
;.re  initial  evidence  for  the  coda  3-factor  correction  of 
equation  (III-17)  was  obtained  from  an  attemot  to  calculate 


consistent  magnitude  estimates  from  measurements  of  frequency- 
wave  number  peaks . 

An  example  is  shown  on  Table  1 1 1 - 1 ,  where  data  were  ex¬ 
tracted  from  measurements  of  spectral  wave  number  peaks  de¬ 
tected  by  the  program  FKCOMB  by  Smart  (1972).  There,  the 
peaks  were  measured  with  varying  time  windows  of  3.2,  6.4, 
and  12. S  seconds.  The  data  were  treated  as  a  stationary 
random  process  with  each  spectral  component  averaged  over  the 
respective  time  window.  The  effect  observed  was  consistent 
fading  of  the  average  amplitude  as  progressively  longer  time- 
windows  were  analyzed.  Note  on  the  bottom  of  Table  1 1 1 - 1 
that  the  B-factor  increases  by  0.1  for  each  doubling  of  the 
time  window  used  to  sample  the  earthquake. 

Consistent  with  observing  a  number  of  events,  the  re¬ 
lationship  of  equation  (III-17)  is  a  simple  coda  decay  model 
reflecting  these  results.  The  parameters  a  and  c  in  equation 
(III-17)  were  estimated  at  a  =  0.33  and  c  =  0.15.  This  im¬ 
plies  a  model  where  the  coda  amplitude,  A  (t)  drops  to  about 
0.7  of  the  signal  amplitude  after  duration,  D,  and  decays  with 
the  dimensionless  factor  (t/D)  raised  to  the  -0.33  power. 

This  is  shown  by  equation  (III-13)  . 

Ac(t)  =0.7  (t/D)'0-'3,  (t >D)  .  (Ill- IS) 

This  decay  model  was  observed  to  be  consistent  with  modified 
3-factors  computed  empirically  for  determining  magnitudes 
from  spectral  wave  number  peaks.  It  is  considered  applicable 
for  the  decay  immediately  following  a  signal  of  duration,  D. 
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TIME  WINDOW  EVALUATION 
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rurtr.er  observed  with  FKCOMB  that  spectral  wavenumbers 
trie  coda  rollowing  a  signal  peak  indicated  frequency  peaks, 
muons ,  and  aT/dds  consistent  with  those  observed  at  the 
oial  signal  peak.  This  provides  some  support  for  the  shot 
se  model  to  describe  the  initial  strong  motion  character- 
ic  of  earthquake  coda. 


Onset  Time  and  Duration  of  a  Propagated  Phase 


Part  of  the  random  characteristics  of  signals  is  taken 
care  or  by  the  assumption  of  the  sudden  onset  of  clustered 
signal  peaks  persisting  for  duration  D,  followed  by  a  decay- 
~?*g  cocia’*  Pot  more  completeness  we  need  to  also  consider  the 
emergence  time  or  a  signal  onset.  Since  we  are  detecting 
signals  on  segmented  records,  there  is  no  need  to  continuous  1, 
update  the  noise  statistics.  Therefore,  there  is  no  problem 
0”  missing  emergent  events \  only  a  problem  of  accurately 
timing  the  beginning  of  such  events.  A  criteria  for  the 
model  or  an  emergent  signal  is  a  sequence  of  monotonically 
-“C*  eas  mg  envelope  peaks.  Also,  a  model  for  the  end  of  a 
signal  duration  D  is  a  sequence  of  monotonically  decreasing 
signal  peaks.  Thus,  at  the  onset  of  a  possible  signal,  the 
start  or  the  time  window  is  held  stationary  as  monotonically 
ascending  peaks  are  sensed.  If  the  maximum  peak  of  such  a 
sequence  passes  signal  detection  threshold  criteria,  then  the 
starting  point  of  the  emergent  signal  is  correctly  ascertained. 
Conversely,  the  end  of  the  signal  window  is  similarly  detei- 
mmec  by  similarly  operating  on  a  monotonically  decreasing 
sequence  or  envelope  peaks.  If  the  minimum  envelope  peak  of 
sucn  a  sequence  passes  signal  rejection  criteria,  then  the 
end  or  the  signal  window  is  accurately  ascertained.  In  short, 


Hi 
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->'  treating  monotonic  envelope  sequences  as  single  maximum 
or  minimum  envelope  peaks,  signal  windows,  including  those 
of  emergent  events,  should  be  accurately  timed. 

5.  Measurement  of  Propagated  Phases 

After  detecting  a  propagated  phase  as  a  cluster  of 
envelope  peaks,  we  determine  a  time  window  of  onset  time,  t, 
ar.c  duration  D.  The  detection  status  and  the  magnitude,  fre¬ 
quency  and  bandwidth  of  each  spectral  moment  is  determined 
rcr  each  detected  signal.  Network  processing  consists  of 
retrieval  of  this  information  for  all  detected  phases  on  the 
records  associated  with  the  preliminary  event  location  under 
consideration.  Absorption  corrections  can  be  performed  at 
some  of  those  stations  where  signals  are  detected  at  anomal¬ 
ously  low  frequencies.  Unless  these  absorption  corrections 
are  applied  and  are  sufficient  to  obtain  frequency  band  con¬ 
vergence,  the  use  of  uncorrected  data  could  reduce  the 
earthquake/explosion  discrimination  power  of  the  data.  Next, 
time  anomaly  and  magnitude  consistency  criteria  can  be  applied 
to  retme  the  process  of  associating  propagating  phases  and 
updating  the  location  of  the  event.  Finally,  the  measure¬ 
ments  cf  signals  and  noise  can  be  analyzed  to  determine  un¬ 
biased  estimates  of  the  magnitude  of  the  event  and  event  dis¬ 
criminants  . 

G.  NETWORK  STRATEGY  FOR  RETRIVEING  SIGNALS  OF  SMALL  EVENTS 

e  expect  to  be  taced  with  the  problem  of  identifying 
small  events  which  are  masked  by  seismic  noise.  By  applying 
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__xed  constant  ra_se  a^arm  rate  thresholds  to  constrain  false 
alarms,  we  may  keep  such  false  alarms  to  an  acceptably  low 
ievei»  cut  -*ae  same  time  make  it  nearly  impossible  to  ex¬ 
tract  weak  signals.  In  the  case  of  low  magnitude  signals, 
we  can  employ  a  station  Variable  False  Alarm  Rate  (YVAR) 
strategy  which  raises  the  expected  ratio  of  retrieved  signals 
to  retrieved  false  alarms  to  an  acceptable  level.  Oi  the  pre¬ 
sumption  that  an  identification  is  required  for  all  events 
analyted,  it  is  better  to  have  at  least  one  or  a  few  signals 
combined  with  some  false  alarms  than  to  have  only  noise  mea¬ 
surements  to  identify  the  signal.  This  strategy  is  made 
somewhat  more  sensible  by  the  fact  that  network  validation 
rocedures  will  be  able  to  eliminate  at  least  some  of  the 
alse  alarms  and  that  network  magnitude  determinations  will 
compensate  to  some  extent  for  noisy  determinations  of  the 
event  parameters. 


This  type  of  post -automat ic  association  network  level 
signal  analysis  effectively  reduces  independent  signal  mea¬ 
surements  to  more  accurate  event  measurements.  It  provides 
a  network  view  to  check  the  consistency  of  associated  signals, 
-ign&ls  with  reasonable  identifiable  source  and  propagation 
chi. ac er is t ics  can  be  passed  on  while  other  'false'  signal 
measurements  can  be  eliminated.  Array  and  polarization  mea¬ 
surements  as  well  as  travel  times  can  be  used  in  this  context 
and  weed  out  obvious  local  events  and  to  identify 
core  phases  outside  the  normal  teleseismic  range,  later  tele- 
seismic  phases,  and  regional  phases. 
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m  this  phase  of  applying  network  strategies,  we  define 
our  ODjective  as  that  of  detecting  snail  signals;  possibly 
e\ en  oe_ow  the  most  frequently  occurring  noise  level  or  by 
a?P ->'-ng  mul t ivariate  waveform  recognition  procedures.  It 
-5  expected  that  the  normally  high  automatic  detector  false 
a.arm  rate  will  be  reduced  by  such  multivariate  analysis  of 
tnese  small  signals.  Furthermore,  network  strategies  applied 
"  L*lese  measurements  such  as  measurement  of  travel  time 
anomalies  and  magnitude  deviations  will  further  reduce  such 
*aije  alarms  to  a  level  which  will  allow  accurate  phase 
association  and  location.  As  previously  pointed  out,  VFAR 
threshold  strategies  can  further  reduce  false  alarms  and  im¬ 
prove  network  multivariate  magnitude  measurements  for  event 
discrimination. 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 


A  concise  description  was  given  for 
multivariate  analytic  detector  to  be  used 
short-period  seismic  signals.  The  output 
matic  seismic  signal  editor  could  provide 
to  improve  the  reliability  and  efficiency 
fication. 


the  design  of  a 
as  an  editor  of 
of  such  an  auto- 
additional  data 
of  event  identi- 


One  of  the  functions  performed  by  the  automatic  signal 
editor  is  to  extract  all  of  the  waveforms  which  can  be  pos¬ 
sibly  interpreted  as  signals  from  the  event.  These  extrac¬ 
tions  include  long-period  surface  waves,  snort-period  tele- 
seismic  P  waves,  secondary  phases,  and  regional  phases. 

This  information  could  be  automatically  inserted  into  signal 
measurement,  files  and  accessed  for  additional  interactive 
seismic  processing  by  seismic  analysts  . 

Another  function  of  the  automatic  signal  editor  is  to 
reduce  basic  ground  motion  measuremen.s  to  multivariate  source 
discriminants  such  as  event  ground  motion  or  spectral  moment 
magnitudes,  including  their  dominant  frequency,  bandwidth 
and  complexity.  These  signal  measurements  can  be  clustered 
to  identify  anomalous  events  and  to  calibrate  normal  earth¬ 
quakes  from  various  source  regions.  These  spectral  moment 
magnitudes  can  be  interpreted  physically,  after  appropriate 
attenuation  corrections,  as  conventional  magnitude  measure¬ 
ments  [log  A/ (Dominant  Period))  of  ground  motion;  i.e.,  dis¬ 
placement  potential,  displacement,  velocity,  acceleration  and 
jerk.  These  are  band-limited  whitened  source  data  in  appro¬ 
priate  roll-off  portions  of  the  signal  spectrum.  Their  use 
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s  many  problems,  including  source  sice-scaling  associated 
the  arbitrary  application  of  narrow  bandpass  filters  to 


compute  source  discriminants. 
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